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PREFACE 



It has been almost thirty years since a book was published that was entirely dedi- 
cated to the theory, description, characterization and measurement of the thermal 
conductivity of solids. For example, the excellent texts by authors such as Berman 1 , 
Tye 2 and Carlslaw & Jaeger 3 remain as the standards in the field of thermal con- 
ductivity. Tremendous efforts were expended in the late 1 950’s and 1960’s in rela- 
tion to the measurement and characterization of the thermal conductivity of solid- 
state materials. These efforts were made by a generation of scientists, who for the 
most part are no longer active, and this expertise would be lost to us unless we are 
aware of the great strides they made during their time. 

I became interested in the field of thermal conductivity in the mid 1990’s in relation 
to my own work on the next generation thermoelectric materials, of which the 
measurement and understanding of thermal conductivity is an integral part. 4 In 
my search for information, I found that most of the books on the subject of 
thermal conductivity were out of print. Much of the theory previously formulated 
by researchers such as Klemens 5 and Slack 6 contain considerable theoretical insight 
into understanding thermal conductivity in solids. However, the discovery of new 
materials over the last several years which possess more complicated crystal struc- 
tures and thus more complicated phonon scattering mechanisms have brought new 
challenges to the theory and experimental understanding of these new materials. 
These include: cage structure materials such as skutterudites and clathrates, metal- 
lic glasses, quasicrystals as well as many of new nano-materials which exist today. 
In addition the development of new materials in the form of thin film and super- 
lattice structures are of great theoretical and technological interest. Subsequently, 
new measurement techniques (such as the 3 -u> technique) and analytical models to 
characterize the thermal conductivity in these novel structures were developed. 
Thus, with the development of many new and novel solid materials and new mea- 
surement techniques, it appeared to be time to produce a more current and readily 
available reference book on the subject of thermal conductivity. Hopefully, this 
book, Thermal Conductivity-2004 : Theory, Properties and Applications, will serve 
not only as a testament to those researchers of past generations whose great care in 
experimental design and thought still stands today but it will also describe many of 
the new developments over the last several years. In addition, this book will serve 
as an extensive resource to the next generation researchers in the field of thermal 
conductivity. 

First and foremost, I express great thanks to the authors who contributed to this 
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PREFACE 



book for their hard work and dedication in producing such an excellent collection 
of chapters. They were very responsive to the many deadlines and requirements and 
they were a great group of people to work with. 1 want to personally acknowledge 
my many conversations with Glen Slack, Julian Goldsmid, George Nolas and 
Ctirad Uher, as well as many other colleagues in the field, as I grasped for the 
knowledge necessary to personally advance in this field of research. I also want to 
acknowledge the support and encouragement of my own institution, Clemson Uni- 
versity (especially my Chairman: Peter Barnes), during this editorial and manu- 
script preparation process. 1 am truly indebted to all my graduate students, for their 
contributions to these volumes, their hard work, for the patience and understanding 
they exemplified during the editorial and writing process. A special thanks goes to 
the publishers and editors at Kluwar Press for their encouragement and patience in 
all stages of development of these manuscripts for publication. 1 am indebted to my 
assistant at Clemson University, Lori McGowan, whose attention to detail and 
hard work (copying, reading, filing, corresponding with authors, etc.) helped 
make this book possible. I especially wish to acknowledge my wonderful wife, 
Penny, and my great kids for their patience and understanding for the many hours 
that I spent on this work. 



Terry M. Tritt 
September 18, 2003 
Clemson University 
Clemson, SC 
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Chapter 1.1 



THEORY OF THERMAL CONDUCTIVITY 



Jihui Yang 

Materials and Processes Laboratory, GM R&D Center, Warren, MI, USA 



1. INTRODUCTION 

Heat energy can be transmitted through solids via electrical carriers (electrons or 
holes), lattice waves (phonons), electromagnetic waves, spin waves, or other excita- 
tions. In metals electrical carriers carry the majority of the heat, while in insulators 
lattice waves are the dominant heat transporter. Normally, the total thermal con- 
ductivity k can be written as a sum of all the components representing various 
excitations : 

( 1 ) 

a 

where a denotes an excitation. The thermal conductivities of solids vary dramati- 
cally both in magnitude and temperature dependence from one material to another. 
This is caused by differences in sample sizes for single crystals or grain sizes for 
polycrystalline samples, lattice defects or imperfections, dislocations, anharmonicity 
of the lattice forces, carrier concentrations, interactions between the carriers and the 
lattice waves, interactions between magnetic ions and the lattice waves, etc. The 
great variety of processes makes the thermal conductivity an interesting area of 
study both experimentally and theoretically. 

Historically thermal conductivity measurement was used as a powerful tool for 
investigating lattice defects or imperfections in solids. In addition to opportunities 
for investigating exciting, intriguing physical phenomena, thermal conductivity 
study is also of great technological interest. Materials with both very high and 
very low thermal conductivities are technologically important. High-thermal-con- 
ductivity materials like diamond or silicon have been extensively studied in part 
because of their potential applications in thermal management of electronics. 1-8 
Low-thermal-conductivity materials like skutterudites, clathrates, half-heuslers, 
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chalcogenides, and novel oxides are the focus of the recent quest for high-efficiency 
thermoelectric materials. 9-34 

The aim of this chapter is to review and explain the main mechanisms and 
models that govern heat conduction in solids. More detailed theoretical treatments 
can be found elsewhere. 35-38 Furthermore, only heat conduction by carriers (elec- 
trons) and lattice waves (phonons) at low temperatures will be discussed. It is in 
this temperature region that most of the theoretical models can be compared 
against experimental results. 



2. SIMPLE KINETIC THEORY 



Thermal conductivity is defined as 



K = — 



Q 

vt’ 



( 2 ) 



where Q is the heat flow rate (or heat flux) vector across a unit cross section 
perpendicular to Q and T is the absolute temperature. For the kinetic formulation 
of thermal conduction in gases, let us assume that c is the heat capacity of each 
particle and n is the concentration of the particles. In the presence of a temperature 
gradient VT, for a particle to travel with velocity v its energy must change at a rate 
of 



dE 

~dt 



cv • VT. 



( 3 ) 



The average distance a particle travels before being scattered is vt, where r is the 
relaxation time. The average total heat flow rate per unit area summing over all 
particles is therefore 

Q = — ncr(v ■ v)VT = — —ncrxiwT . (4) 



The brackets in Eq. (4) represent an average over all particles. Combining Eqs. (2) 
and (4), we have 



1 2 V , 

k = -UCTV = -0 vl, 

3 3 



( 5 ) 



where C = nc is the total heat capacity and l = vt is the particle mean free path. In 
solids the same derivation can be made for various excitations (electrons, phonons, 
photons, etc.). Equation (5) can then be generalized to 



K = 



1 

3 






C n 



C : 



( 6 ) 



where the summation is over all excitations, denoted by a. In general, Eq. (6) gives 
a good phenomenological description of the thermal conductivity, and it is practi- 
cally very useful for order of magnitude estimates. 

Like most of the nonequilibrium transport parameters, thermal conductivity 
cannot be solved exactly. Calculations are usually based on a combination of per- 
turbation theory and the Boltzmann equation, which are the bases for analyzing the 
microscopic processes that govern the heat conduction by carriers and lattice waves. 
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3. ELECTRONIC THERMAL CONDUCTION 



The free electron theory of electron conduction in solids in the first instance con- 
siders each electron as moving in a periodic potential produced by the ions and 
other electrons without disturbance, and then regards the deviation from the per- 
iodicity due to the vibrations of the lattice as a perturbation. The possible values of 
the electron wave vector k depend on the periodicity and the size of the crystal. The 
fc-space is separated into Brillouin zones. The electron energy Eg depends on the 
form of the potential and is a continuous function of k in each zone, but it is 
discontinuous at the zone boundaries. The values of Eg in a zone trace out a 
"band” of energy values. 

The distribution function that measures the number of electrons in the state k 
and at the location r is fg. In equilibrium, the distribution function is /'!, given by 




where E F and are the Fermi energy and Boltzmann's constant, respectively. 
According to the Boltzmann equation, in the presence of an electrical field E 
and a temperature gradient VT, the steady state that represents a balance between 
the effects of the scattering processes and the external field and temperature gra- 
dient can be described as 



fjzA 

r(fc) 



-Vr 



( d /“ df, f 

k -S7T + e—^E 



dT 



dE t 



( 8 ) 



where r(fc), v g, and e are the relaxation time, the electron velocity, and the electron 
charge, respectively. Taking into account the spatial variation of the Fermi energy 
and the explicit expression of f~, Eq. (8) can be written in the form 



h-fl 

r(fc) 



-Vr 



Eg-E F df, df, 

^ 5 i VT+e si 




( 9 ) 



The effective field acting on the electrons is E e $ - E - w; ' . Equation (9) can be 
entered into the expressions for the electrical current density J and the flux of 
energy Q : 

J = j evgfgdk (10) 



and 



Q = 



(Eg-E F )vgfgdk. 



( 11 ) 



If we define a general integral K n as 



I<n 




e vg) 2 r(k)(Eg-E F ) n 



d A 

dEi 



dk. 



( 12 ) 



then Eqs. (10) and (11) can be written as 
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and 



J = e 1 K 0 E e f f — — K\WT 


(13) 


Q = eKiE e ff - l- K 2 WT , 


(14) 



respectively. The electronic thermal conductivity can be found with J = 0 such that 



He = - 



Q 



VT 



1 



j— o 



T 



= - K 2 --^ . 



K\ 



K c 



(15) 



(16) 



The electrical conductivity a can be derived from Eq. (13) as 

a = e 2 Ko . 

Since QE_ is approximately a delta function at the Fermi surface with width 
ksT, A"„/can fee evaluated by expansion: 

a{E, 



Kn = 

Therefore 39 

and 



' t7 V?l -) , ^ ,„2 rn2 & 

e 2 



( E u-E F ) n ^+‘-k 2 B T 2 — 2 



{E S -E F )"-^ 



tt k 2 B T 2 

Kl ~ 



<t(E f ) + o(^bTL\ 






Er=E F ■ 



(17) 



(18) 



(19) 



Equation (15) through (19) lead to the Wiedemann-Franz law with the standard 
Lorentz number La as 



K e 7T 2 k? B 

0 _ crT ~ 3 e 2 ■ 



( 20 ) 



The numerical value of Lq is 2.4453 x 10~ 8 W-fl/K 2 . This shows that all metals have 
the same electronic thermal conductivity to electrical conductivity ratio, and this 
ratio is proportional to the absolute temperature. Furthermore, for strong degen- 
erate electron gases, Eq. (20) is independent of the scattering mechanism and the 
band structure for the electrons as long as the scatterings are elastic. The Wiede- 
mann-Franz law is generally well obeyed at high-temperatures. In the low- and 
intermediate-temperature regions, however, the law fails due to the inelastic scat- 
tering of the charge carriers. 36, 37, 39 

Typically one needs to calculate the relaxation times for the various relevant 
electron scattering processes and use Eqs. (15) and (16) to determine the electronic 
thermal conductivity and the electrical conductivity. Over a wide temperature 
range, scattering of electrons by phonons is a major factor for determining the 
electrical and electronic thermal conductivities. The resistance due to this type of 
scattering is called the ideal resistance. The ideal electrical and electronic thermal 
resistances can be approximately written as 36 
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Pi = 




and 



W t = 



1 



Ki 



A 

ut 





1 Ji(O d /T)\ 
2n 2 Js{9 D /T)y 



respectively, where 





(e x - l) 2 



dx , 



( 21 ) 

( 22 ) 

(23) 



3t rhg^(G') 2 
4 e 2 (m*) 2 n c k B 0 D k 2 F v 2 F ’ 



(24) 



On is the Debye temperature for phonons (to be discussed below), hr is the electron 
wave number at the Fermi surface, qu is the phonon Debye wave number, h is the 
Planck constant, G' is the constant representing the strength of the electron-pho- 
non interaction, m* is the electron effective mass, n c is the number of unit cells per 
unit volume, and vf is the electron velocity at the Fermi surface. According to Eqs. 
(21) through (24), at high-temperatures (T >> 0 D ) 



Pi = 



AT 

49 ^ 



(25) 



and 



Wi = 



A 

49fLo 



Pi 

L 0 T - 



(26) 



The Wiedemann-Franz law is obeyed; W t is a constant. At low temperatures 
(T« 0 D ) 

Pi = 124. 4A(-^-) 5 (27) 



and 



Wi = 124.4 



A /T\ 3 3 



L 0 T \e D J ir 2 




(28) 



In addition to the electron-phonon interaction, the electron-defect interaction con- 
tributes to the electrical resistivity and electronic thermal resistivity. One can ap- 
proximately write the resistivity p and electronic thermal resistivity W e as (Mat- 
thiessen’s Rule) 

p = -= po + Pi (29) 

a 



and 

W e = — =W 0 + Wi. 

K e 



( 30 ) 
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Here; p 0 and Wo are the residual electrical resistivity and electronic thermal resis- 
tivity, respectively, caused by electron scattering due to impurities and defects, and 
po is independent of temperature. If we assume that all defects scatter electrons 
elastically, p 0 and Wo should be related by the Wiedemann-Franz law: 

Wo = po/LqT ex l/T. (31) 



At very low temperature, since Wo increases and IF, decreases with decreasing 
temperature, 



1 L 0 
Wo po 



(32) 



The electronic thermal conductivity increases linearly with increasing temperature. 
As temperature increases, IF becomes relatively more important. If W, increases to 
values comparable to that of Wo at sufficiently low temperature, then K e will pass a 
maximum, decrease, and eventually attain the high-temperature constant described 
by Eq. (26). In alloys or metals with a high concentration of defects, W 0 and IF, 
only become comparable at high-temperatures. In this case, there is no maximum in 
the K e , versus, T curve, and n, will approach the high-tenrperature constant rnono- 
tonically as T increases. Both cases are illustrated in Fig. 1, which shows the low- 
temperature thermal conductivity of silver measured by White. 40 In Fig. 1, sample 
(a) is obtained by annealing sample (b). Low, temperature thermal conductivity is 
therefore dramatically affected by the residual resistance. By combining Eqs. (27) 
through (31), we have 




0 20 40 60 80 100 120 140 



T (K) 

FIGURE 1 Thermal conductivity of two silver samples with (a) very low and (b) moderately high 
residual resistivities from Ref. 40, showing the profound influence of the residual resistivity on the 
low-temperature thermal conductivity. Sample (a) was obtained by annealing sample (b). Reprinted 
with permission from IOP Publishing Limited. 
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T/e D 

FIGURE 2 Calculated (K e /(7T)/I,o versus T/9d curves for various values of pal A, and for monovalent 
metals. The calculations are based on Eq. (33) using (kp/qp,) 2 = ^/[3]2/2 for monovalent metals. 



f+0 5 ^) 



e /<jT _ 

Lo p a + © 5 {[1 + ^(|) -2 



(33) 



This quantity versus T/0 D is shown in Fig. 2 for various values of and for 
monovalent metals ((/cf/?d) 2 = \/2/2). Figure 2 shows that the Wiedemann-Franz 
law holds at very high and very low temperatures except for very pure metals, but 
K e /aT underestimates L 0 at intermediate temperatures in a manner that strongly 
depends on the amount of impurity present. 

It should be stressed that Eqs. (21) and (22) were derived without electron- 
phonon umklapp processes, and the strength of the electron-phonon interaction 
was assumed to be a constant. Close agreement between the model and experimen- 
tal data should not be expected in most cases. Modifications to the model are 
discussed elsewhere. 37 

The discussion so far has been focused on degenerate electron gases (like those in 
metals). In semiconductors, because of the energy gap Eg between the top of the 
valence band and the bottom of the conduction band and the consequence that the 
Fermi level lies in the gap, the electron distribution may no longer be degenerate. 
Also, there may exist positively charged particles, the “holes,” that may also con- 
tribute to the transport. In this case, the distribution functions for electrons and 
holes may be written as 



fe 



exp 



E e + Ep 
k B T 



(34) 



and 



f° 

Jh 



exp 



E h + Ep 
k B T 



(35) 



where and are the equilibrium distribution functions for electrons and holes, 
respectively, E e and E h are the electron and hole energies, respectively, and Ep and 
Ep are the distances of the Fermi level below the bottom of the conduction band 
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and above the top of the valence band, respectively. If one wishes to use the general 
integral defined in Eq. (12), the conductivities can be written similarly to those in 
Eqs. (15) and (16): 



K e 



1 

T 



(K 2 + 1<' 2 ) 



{Ki + K\f[ 

K 0 + K' 0 j 



(36) 



and 

er= e 2 (K 0 + K'q), 



(37) 



where K n and K' n denote integrals for electrons and holes, respectively, and the 
distribution functions for the integrals should be those listed in Eqs. (34) and (35). 
In general, the ratio K e /aT is rather complicated to work out. For simplicity, 
assume that the energy bands are parabolic and the energy dependences of the 
carriers’ relaxation times are the same such that 

Te,h(E) cx E a , (38) 



where a is a constant, and the subscripts e and h denote the electrons and holes. 
With some algebra one finds 36 



K e 

aT 




^5 + 2a + 



Eg \ 2 n e iJ, e n h iJ, h 
k bT J (n e He + rihUh) 2 



(39) 



In Eq. (39), n and p are the carrier concentration and mobility, and again the 
subscripts e and h denote electrons and holes. The first term in Eq. (39) is the 
standard Lorentz number for nondegenerate Fermi distributions. When carrier- 
acoustic phonon scattering dominates, a = -1/2, and if n e = nu, Eq. (39) can be 
written as 



K e = aT 




2 + 




Eg \ 2 He/Hh 
kBTJ (1+Pe/V) 2 



(40) 



The second term in Eq. (39) or (40) is called the bipolar diffusion term. This type of 
energy transport in addition to that carried by the electrons and the holes is due to 
the creation of electron-hole pairs that extract an amount of energy Eg at the high- 
temperature end; on recombination this energy is given up to the cold end. The 
bipolar diffusion may be a noticeable contribution to the electronic thermal con- 
ductivity when both the carrier concentration and the mobility are about equal for 
electrons and holes, and their mobilities have reasonably high values. 41,42 This 
should be the case for intrinsic semiconductors. In the case of doped semiconduc- 
tors only one of the carriers has high mobility; the low-mobility carrier will not be 
fast enough to accompany the high-mobility carrier for the recombination at the 
cold end; therefore, the bipolar contribution may not be noticeable. Anomalously, 
high thermal conductivity at high-temperatures for InSb, Ge, Si, and Bi may be 
explained by this bipolar diffusion. 43, 44, 45 

Figure 3 shows the temperature dependence of (IT - IT /) /T for pure Ge mea- 
sured by Glassbrenner and Slack. 44 The thermal resistance due to isotopes IT/ is 
almost independent of temperature and is subtracted from the total thermal resis- 
tance W. The term (IT — Wi)/T increases linearly as a function of T for 7' < 700 
K, reaches a maximum at about 700 K, and decreases with increasing T for T > 
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FIGURE 3 The temperature dependence of (W — Wi)/T for pure Ge from Ref. 44. W is the total 
measured thermal resistivity, and W/ is the calculated thermal resistance due to isotopes. The dashed 
line is a theoretical extrapolation for thermal resistivity induced by phonon-phonon interactions. 



700 K. The theoretical extrapolation (dashed line) represents T and T 2 components 
of the thermal resistivity due to phonon-phonon interactions. The former is due to 
three-phonon umklapp processes (to be discussed), and the latter is attributed to 
four-phonon processes. 46 Deviation from the dashed line at high-temperatures is 
ascribed to the electronic thermal resistance. It was found that the contribution to 
the total thermal conductivity from the first term of Eq. (40) is less than 10% up to 
the melting point. The majority of the deviation between (W - Wi)/T and the 
dashed line at high-temperatures is due to dipolar diffusion. The values of Eq 
can be estimated by subtracting the lattice component and the first term of Eq. 
(40) from the total thermal conductivity and comparing it against the second term 
in Eq. (40). The results are in agreement with the published results. 44 

3.1. Lattice Thermal Conductivity 

Lattice thermal conduction is the dominant thermal conduction mechanism in non- 
metals, if not the only one. Even in some semiconductors and alloys, it dominates a 
wide temperature range. In solids atoms vibrate about their equilibrium positions 
(crystal lattice). The vibrations of atoms are not independent of each other, but are 
rather strongly coupled with neighboring atoms. The crystal lattice vibration can be 
characterized by the normal modes, or standing waves. The quanta of the crystal 
vibrational held are referred to as “phonons.” In the presence of a temperature 
gradient, the thermal energy is considered as propagating by means of wave packets 
consisting of various normal modes, or phonons. Derivation of phonon dispersion 
curves (wj versus q curves, where wj and q are the phonon frequency and wave 
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FIGURE 4 Schematic phonon dispersion curves for a given direction of q of (a) monatomic lattice and 
(b) diatomic lattice. The lattice parameter is denoted ao. 



vector, respectively) can be found in a standard solid-state physics book. 47 Sche- 
matic phonon dispersion curves for monatomic and diatomic lattices are shown in 
Fig. 4. Phonon dispersion curves for solids normally consist of acoustic and optical 
branches. The low-frequency acoustic branches correspond to atoms in a unit cell 
moving in same phase, whereas the high-frequency optical branches represent 
atoms in a unit cell moving in opposite phases. Normally optical phonons them- 
selves are not effective in transporting heat energy because of their small group 
velocity dco^/dq, but they may affect the heat conduction by interacting with the 
acoustic phonons that are the main heat conductors. 

The phonon distribution function, which represents the average number of pho- 
nons with wave vector q, is N§. In equilibrium, the phonon distribution function 
can be written as 

/vr2 = 1 141) 

q exp(hujj/kBT) — 1 



The Boltzmann equation assumes that the scattering processes tend to restore a 
phonon distribution to its equilibrium form Ni at a rate proportional to the 
departure of the distribution from equilibrium, such that 



Nq-N l 



— (tig • VT) 



dN% 



(42) 



where v g is the phonon group velocity and r q is the phonon scattering relaxation 
time. The heat flux due to a phonon mode q is the product of the average phonon 
energy and the group velocity. Therefore the total heat flux carried by all phonon 
modes can be written as 

Q = ^2 NqhiOqVg. (43) 

9 



Substituting Eq. (42) into Eq. (43) yields 

dNt 



Q = — (cos 2 9) ■ 



dT 



i am. 

VT = --£^7^VT s 



3-V 9 g q dT 

Q 



( 44 ) 



where 9 is the angle between v g and VT; the lattice thermal conductivity is 
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KL = - 



Q_ 

VT 



« 



T Q 



dNi 

~df ' 



(45) 



At this point approximations need to be made to Eq. (45) in order to obtain 
meaningful results. Furthermore, evaluations of various phonon scattering relaxa- 
tion times are usually difficult to make precise. It is not worthwhile to try to 
calculate Eq. (45) for the precise phonon frequency spectrum and dispersion curve 
of a real solid. Assumptions of Debye theory should be used: an average phonon 
velocity v (approximately equal to the velocity of sound in solids) is used to replace 
v g , ujg = vg for all the phonon branches, and the phonon velocities are the same for 
all polarizations. The summation in Eq. (45) can be replaced by the integral 







cW2 

f(q)dq, 



(46) 



where f(q)dq = (3g 2 /27r 2 ) dq, and therefore f(u)cL> = (3cu 2 /27r 2 t> 3 ) dw. Using the 
Debye assumptions and Eqs. (41) and (45) leads to 



kl 




flUT T q (u>) 



{Too / k B T 2 ) exp(Tiu> / k B T) 

— dw. 

[exp(hu> / k B T) — 1]" 



(47) 



where u B is the Debye frequency such that 



3 N = 



,f(uj)du 



(48) 



is the total number of all distinguishable phonon modes. If we make the substitu- 
tion x = Tiuj/k B T and define the Debye temperature d D = huj D /k B , Eq. (47) be- 
comes 



kl 



Kb ,fa'3 T 3 
ln 2 v h 




x 4 e x 

(■ e x — l) 2 



;dx. 



(49) 



Within the Debye approximation the differential lattice specific heat is 



x 4 e x 






-dx. 



(50) 



If we define the mean free path of the phonons as l(x) = vr q (x), the lattice thermal 
conductivity can be written as 

i rOo/T i rOo/T 

kl = - / v 2 T q {x)C{x)dx = - / C(x)vl(x)dx. (51) 

3 Jo 3 Jo 



This is analogous to the thermal conductivity formula [Eq. (6)] derived from simple 
kinetic theory. 

Equation (49) is usually called the Debye approximation for the lattice thermal 
conductivity. If one can calculate the relaxation times t,(x) for various phonon 
scattering processes in solids and add the scattering rates such that 

Tg\x) = ^rr’Or), 



(52) 
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Eq. (49) should be sufficient to obtain the lattice thermal conductivity. It is indeed 
then adequate for analyzing and predicting a lot of the experimental data especially 
when a large concentration of lattice defects prevails in the solid. Sometimes, than 
even in the case of some chemically pure crystals, because of the existence of 
appreciable amounts of isotopes, the lattice thermal conductivity can be represented 
satisfactorily by the Debye approximation. 

The phonon scattering processes included in the Debye approximate are resistive 
and are called umklapp processes or U-processes. The total crystal momentum is 
not conserved for U-processes. Because such processes tend to restore non-equili- 
brium phonon distribution to the equilibrium distribution described by Eq. (42), 
they give rise to thermal resistance. There exist, however, other nonresistive and 
total crystal-momentum-conserving processes that do not contribute to the thermal 
resistance but may still have profound influence on the lattice thermal conductivity 
of solids. Such processes are called normal processes or N-processes. Even though 
N-processes themselves do not contribute to thermal resistance directly, they have 
the great effect of transferring energy between different phonon modes, thus pre- 
venting large deviations from the equilibrium distribution. 

Since the N-processes themselves do not tend to restore the phonon equilibrium 
distribution, they cannot be simply added to Eq. (52). The Callaway model is the 
most widely used in analyzing the effect of N-processes on the lattice thermal 
conductivity. 48 Callaway’s model assumes that N-processes tend to restore a non- 
equilibrium phonon distribution to a displaced phonon distribution of the form 49 

7V-(X) = L = Nl + \ L| «p(^A»r) (53) 

exp[(/ku —5 -A )/k B T] — 1 k B T [exjp(hLd/k B T) — 1] 

where A is some constant vector (in the direction of the temperature gradient) that 
determines the anisotropy of the phonon distribution and the total phonon mo- 
mentum. If the relaxation time for the N-processes is t n , Eq. (42) can be modified 
to 



| Ng-Ng( A) 

Tq Tn 



-(v-VT) 



dN9 

~dT' 



(54) 



If we define a combined relaxation time r c by 



T c 1 — T q * + t n > 




the Boltzmann equation [Eq. (54)] can be written as 

_ huj ^ exp (hui/k B T) q- A exp(Hcd/k B T) _ n\ 

k B T 2 [exp(/z u>/k B T) — 1]" TNk B T [exp (hui/k B T) — 1]~ t c 



If we express n\ as 



v huJ gt cxp(Mar) 

q k B T 2[ [exp^hcd/ksT) — l ] 2 ' 



then Eq. (57) can be simplified to 
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fi "“« £ :.vr+i^ = ^s.VT. 



T r T 



tn 



T 



(59) 



Since A is in the direction of the temperature gradient, it is convenient to define a 
parameter j3 with the same dimension as the relaxation time: 

A=-J/3 v 2 VT. (60) 



T 



Because q = vuj/v 1 , Eq. (59) can be further simplified to 

Oiq = r c (l + 0/t n ). 



(61) 



From Eq. (58) it is straightforward [the same procedure as Eqs. (42) through (49)] 
to show that the lattice thermal conductivity can be expressed as 



k B ( k B ,3rp3 



i-Sd/T 



l( x ) 



x 4 e x 



( e x — l ) 2 



; dx 



2ir 2 v 






c (a:)(l 



0 



x 4 e x 



tn(x) (e 1 — 1)' 



dx. 



(62) 



The task now is to determine (3. Because the total crystal momentum is conserved 
for the N-processes, the rate of phonon momentum change is zero. Therefore 



■ N$-Nn 

tn 



qdq = 0. 



(63) 



Substituting Eqs. (53) and (58) into Eq. (63), we have 
f exp(/fcj fk B T) 



[exp(/i co/k B T) — 1]“ k B T 2 



2 a q {v-S7T) + f-^] dq=0 . 



k B T t n 



(64) 



This can be further simplified by using Eqs. (60) and (61), so that 

[ exp (fiuj/k B T) Huj .. vu> 

1 v ' ’ -{v-VT){a q - /3)—^dq = 0. 



[exp(/z u>/k B T) — l] 7 k B T 2 



T N v- 



(65) 



By inserting Eq. (63) into Eq. (65) and using the dimensionless x as defined earlier, 
we can solve for /?: 



0 = 



rSD/T T r e X X 4 



r o Tv (e x - 1) 



=rdx 



rSn/T 



e x x 4 



1 0 T NTq (e x — 1) 



=rdx. 



( 66 ) 



Therefore the lattice thermal conductivity can be written as 

KL = Kl + K2, 



(67) 



where 



« l = 



k B f k B 

2n 2 v \ h 



r^o/T 



T i 



x 4 e x , 

t c ~dx 

(e x - l) 2 



(67a) 



and 
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K2 



k B 

lir 2 v 




3 

T 3 



/ rSo/T Tc x 4 e 1 
Uo t n ( e *-l ) 2 

rfo/T Tc X i e x 
JO rjvr, ( e i_l) 



dx)~ 



jdx 



(67 b) 



When the impurity level is significant and all phonon modes are strongly scattered 
by the resistive processes in a solid, then t n » r q and t c ~ r 7 . Under this circum- 
stance n\ » k 2 and k l is given by Eq. (67a), the same as Eq. (49) since the re- 
processes do not exist. In the opposite extreme, when N-processes are the only 
phonon scattering processes, we have tn «r q and r c « tn- The denominator of 
K 2 then approaches 0, leading to infinite lattice thermal conductivity as expected 
because the N-processes do not give rise to thermal resistance. 

Now that we have derived the formula for the lattice thermal conductivity, the 
problem is to calculate the relaxation times. Phonon scatterings have been treated 
by numerous authors. Here we list the main conclusions. 

For phonon-phonon normal scattering the relaxation rate 

r^ 1 = Bu> a T b (68) 



is the general form suggested by best fits to experimental thermal conductivity data; 
B is a constant independent of uj and T, ( a , b) - (1,3) was recommended for LiF 
and diamond 50 ’ 51 , and (a, b) = (1, 4) and (2, 3) were used for some group IV and 
III-V semiconductors. 8 
Peierls suggested the form 

Tj) 1 cx T n exp(0£)/mT) (69) 



for the phonon-phonon umklapp scattering with constants n and m on the order of 
l. 52 Based on the Leibfried and Schlomann model, 53 Slack et al. proposed the 
following form for the Griineisen constant 7 and the average atomic mass of M 
in the crystal: 54 

fory 2 

T v ™M^^ Texp (- 0D/3T) - (70) 



Other empirical n and in values were also used, 50,55-57 all of which were based again 
on best fits to experimental data. At sufficiently high-temperatures kl(xI/T if 
phonon-phonon umklapp scattering is the dominant process. 

Klemens was the first to calculate the relaxation rate for phonon-point-defect 
scattering where the linear dimensions of the defects are much smaller than the 
phonon wavelength. 58 The corresponding phonon-point-defect scattering rate is 



t pd 



V 

47 rv 3 



l E/.( 



m — m 



i\2 



m 



(71) 



where V is the volume per atom, m, is the mass of an atom, /, is the fraction of 
atoms with mass m it and rh is the average mass of all atoms. A strain field mod- 
ification to Eq. (71) has been described by Abeles. 59 

The phonon-boundary scattering rate is independent of phonon frequency and 
temperature and can be written as 

r -1 = v/d, 
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where d is the sample size for a single crystal or the grain size for a polycrystalline 
sample. 

For phonon-dislocation scattering, Nabarro separated the effects of the core 
from the surrounding strain field. 60 The corresponding relaxation rates are 



t cL (xNd v 2^ 



(72) 



and 



T Str « N D 



rBjju 
2tt ’ 



(73) 



where N D is the number of dislocation lines per unit area, r is the core radius, and 
B D i is the Burgers vector of the dislocation. 

Pohl suggested an empirical nonmagnetic phonon-resonance scattering relaxation 
rate of 



^ Res 



Cco 2 



(74) 



where C is a constant proportional to the concentration of the resonant defects and 
wo is the resonance frequency. 61 This formula accounted well for the observed low- 
temperature dip in the thermal conductivity of KNCb-containing KC1 crystals. It 
has also been used for fitting experimental data for clathrates and skutterudites. 22,23 

According to Ziman, the relaxation time for the scattering of phonons by elec- 
trons in the conduction state is given by 62 

_i _ e 2 {m*fv ( k B T ^ 

T E pc 4 n % 4 d \\m*v 2 J 



hw 1 + exp[(Im*u 2 - E F )/k B T + h 2 w 2 /8m*u 2 k B T + hw/2k B T]l 

k B T 1 + exp[(^m*w 2 — Ep) /k B T + h 2 w 2 /8m*u 2 k B T — hw/2k B T] J ’ 



where e is the electron-phonon interaction constant or deformation potential and d 
is the mass density. The relaxation time for scattering of phonons by electrons in a 
bound state, as given by Griffin and Carruthers, 63 is 

_j Gw 4 1 

TEPB ~~ [w 2 - (4A/?i) 2 ] 2 [1 + r 2 cu 2 / 4u 2 ] 8 ’ 

where G is a proportionality constant containing the number of scattering centers, 
A is the chemical shift related to the splitting of electronic states, and r 0 is the mean 
radius of the localized state. It should be pointed out that these formulas for 
electron-phonon interaction are based on the adiabatic principle and perturbation 
theory. As argued by Ziman, 36 the theory is only valid if the mean free path of 
electrons L e satisfies the condition 

1 <qL e . (77) 

Since the phonon wave vector q = 2n/X p , this means that the wavelength of a 
phonon X p must not be longer than the mean free path of the electron it scatters. 
The model originally developed by Pippard for explaining the ultrasonic attenua- 
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tion in metals is, however, applicable over the entire range of qL e . M Pippard's 
relaxation times are 



t ep 



4nm*VFL e ur 
15 dv 2 



for qL e « 1 



(78) 



and 



t ep ~ 



7 jnm*VFio 

6 dv 



for qL e >> 1 



(79) 



where n is the electron concentration and vf is the Fermi velocity. There have been 
several cases in which Pippard’s theory was applied to lattice thermal con- 
ductivity. 21,65,66 

Typically it is necessary to use the full Callaway model [Eqs. (67), (67a), and 
(67b)] to interpret experimental data on the effect of isotopes when one starts with 
an isotopically pure crystal. The addition of a small amount of defects will rapidly 
suppress K 2 , which becomes negligible for an impure sample. The decrease of k\ 
upon increasing defect concentration is much slower. For samples with an appreci- 
able amount of defects (even isotopic defects) it is sufficient to use the Debye 
approximation [Eq. (49)] for examining the low-temperature lattice thermal con- 
ductivity. The measured thermal conductivity versus temperature curve is usually fit 
by trial and error for one sample, where an appropriate relaxation rate is carefully 
chosen for each scattering mechanism believed to be present. The same curves for 
samples with additional defects are then fit by using suitable relaxation rates to 
reflect scattering by additional defects. Much interesting physics has been revealed 
by this method. Different phonon scattering processes usually dominate in different 
temperature ranges. Figure 5 plots the lattice thermal conductivity of a binary 




FIGURE 5 The lattice thermal conductivity versus temperature of a CoSb 3 sample from Ref. 21. The 
dots and the solid line represent the experimental data and the theoretical fit. The dashed curves are the 
theoretical limits imposed on the phonon heat transport by boundary scatterings, by a combination of 
boundary plus point-defect scatterings, and by umklapp scatterings. 
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skutterudite compound CoSb 3 . 16 The dots and the solid line represent experimental 
data and a theoretical fit using the Debye approximation, respectively. The possible 
phonon scattering mechanisms are phonon-boundary, phonon-point defect, and 
phonon-phonon umklapp scatterings. The dashed lines in Fig. 5 correspond to 
the theoretical limits on the lattice thermal conductivity set by boundary scattering, 
a combination of boundary plus point-defect scattering, and umklapp scattering. At 
high-temperatures (close to 9d = 300 K), umklapp is the dominant phonon scatter- 
ing mechanism, while boundary and point-defect scattering dominate at low and 
intermediate temperatures, respectively. 

In practice, the dominant phonon method works quite well in predicting the 
effect of phonon scattering processes. The dominant phonon method assumes 
that at a given temperature all phonons are concentrated about a particular domi- 
nant frequency w dom ~ k B T/h. An empirical power law can then be deduced as 
follows. If one particular defect can be described by 1/r oc uf or, equivalently, by 
1/r oc x a T a , taking C oc T 3 (for low T) and employing the simple kinetic formula 
[Eq. (5)] yields kl oc T 3_a . Even though the dominant phonon method is not math- 
ematically justified, the power law is usually valid at low temperatures. For exam- 
ple, one should have n L oc T 3 at low temperatures if boundary scattering is the 
dominant phonon scattering mechanism. 



4. SUMMARY 

In this chapter theoretical treatments of the electronic and the lattice thermal con- 
ductivities at low temperatures have been reviewed. The electronic thermal resis- 
tance for degenerate electron gases is the sum of residual and ideal components. At 
low temperatures W e « AT + B/T 2 , where A and B are constants, while at high- 
temperatures W e approaches a constant. Except for very pure metals, the Wiede- 
mann-Franz law holds well at low and high-temperatures with standard Lorenz 
number Lq. For intermediate temperatures, the Wiedemann-Franz law breaks 
down in a way strongly dependent on the amount of the impurity. For intrinsic 
semiconductors, the bipolar diffusion process enhances the electronic thermal con- 
ductivity. The N- and U-processes are both important for analyzing and predicting 
the lattice thermal conductivity of solids. The lattice thermal conductivity for iso- 
topically pure crystals can be well described by the full Callaway model. When 
impurity concentrations are appreciable, the Debye approximation is adequate for 
modeling experimental data. The different resistive phonon scattering processes in 
the Debye approximation dominate in different temperature ranges. At low tem- 
peratures, the power law predicted by the dominant phonon method is often useful 
in identifying the characteristics of the phonon scattering processes. 
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1. INTRODUCTION 

Metals represent a vast number of materials that have been the backbone of in- 
dustrial development during the past two centuries. The importance of this devel- 
opment upon future technological progress is undiminished and unquestionable. 
Whether in pure, elemental form or as new lightweight, high-strength alloys, metals 
are simply indispensable to modern industrial society. Metals are, of course, known 
for their lustrous and shiny appearance, for their malleability and ductility and, 
above all, for their ability to conduct electric current. Because of their myriad 
applications, it is highly desirable to be able to tailor the properties of metals to 
match and optimize their use for specific tasks. Often among the important criteria 
is how well a given metal or alloy conducts heat. 

The physical parameter that characterizes and quantifies the material’s ability to 
conduct heat is called thermal conductivity, often designated by n. Understanding 
the nature of heat conduction process in metals and being able to predict how well 
a particular alloy will conduct heat are issues of scientific and technological interest. 
In this section we review the fundamental physical principles that underscore the 
phenomenon of heat conduction in metals, develop an understanding of why some 
metals are better than others in their ability to conduct heat, and illustrate the 
behavior of thermal conductivity on several specific examples encompassing pure 
metals and alloys. 
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1.1. Carriers of Heat in Metals 

Metals are solids* and as such they possess crystalline structure where the ions 
(nuclei with their surrounding shells of core electrons) occupy translationally 
equivalent positions in the crystal lattice. We know see Chapter (1.1) that crystal- 
line lattices support heat flow when an external thermal gradient is imposed on the 
structure. Thus, like every other solid material, metals possess a component of heat 
conduction associated with the vibrations of the lattice called lattice (or phonon) 
thermal conductivity, k p . The unique feature of metals as far as their structure is 
concerned is the presence of charge carriers, specifically electrons. These point 
charge entities are responsible not just for the transport of charge (i.e., for the 
electric current) but also for the transport of heat. Their contribution to the thermal 
conductivity is referred to as the electronic thermal conductivity n e . In fact, in pure 
metals such as gold, silver, copper, and aluminum, the heat current associated with 
the flow of electrons by far exceeds a small contribution due to the flow of pho- 
nons, so, for all practical purposes (and essentially for the entire regime of tem- 
peratures from sub-Kelvin to the melting point), the thermal conductivity can be 
taken as that due to the charge carriers. In other metals (e.g., transition metals and 
certainly in alloys) the electronic term is less dominant, and one has to take into 
account the phonon contribution in order to properly assess the heat conducting 
potential of such materials. 

In discussions of the heat transport in metals (and in semiconductors), one makes 
an implicit and essential assumption that the charge carriers and lattice vibrations 
(phonons) are independent entities. They are described by their respective unper- 
turbed wave functions, and any kind of interaction between the charge carriers and 
lattice vibrations enters the theory subsequently in the form of transitions between 
the unperturbed states. This suggests that one can express the overall thermal 
conductivity of metals (and other solids) as consisting of two independent terms, 
the phonon contribution and the electronic contribution: 

K = K e + K p . ( 1 ) 

These two — electrons and phonons — are certainly the main heat carrying entities. 
However, there are other possible excitations in the structure of metals, such as spin 
waves, that may, under certain circumstances, contribute a small additional term to 
the thermal conductivity. For the most part, we shall not consider these small and 
often conjectured contributions. 

It is important to recognize that the theory of heat conduction, whether in 
nonmetallic or metallic systems, represents an exceptionally complex many-body 
quantum statistical problem. As such, it is unreasonable to assume that the theory 
will be able to describe all the nuances in the heat conduction behavior of any given 
material or that it will predict a value for the thermal conductivity that exactly 
matches that obtained from experiment 1 . What one really hopes for is to capture 
the general trend in the behavior of the thermal conductivity, be it among the group 
of materials or in regard to their temperature dependence, and have some reason- 
ably reliable guidelines and perhaps predictive power as to whether a particular 



We consider here metals in their solid form only. 

' Measurements of thermal conductivity are among the most difficult of the transport studies, and the 
accuracy of the data itself is usually no better than 2%. 
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class of solids is likely to be useful in applications where heat carrying ability is an 
important concern or consideration. 

The theory of heat conduction has been outlined in Chapter 1.1. There are 
numerous other monographs and review articles where this topic is treated in detail. 
I am particularly fond of the following texts: Berman’s Thermal Conduction in 
Solids ,' Ashcroft and Mermin’s Solid State Physics , 2 Ziman’s Electrons and 
Phonons , 3 Smith and Jensen’s Transport Phenomena , 4 and Mahan’s Many-Particle 
Physics , 5 The level of coverage varies from an introductory treatment to a com- 
prehensive quantum statistical description. On one hand, a too elementary descrip- 
tion is unlikely to provide much more beyond an outline of the phenomenon; on 
the other hand, a many-body quantum statistical exposure is likely to overwhelm 
most readers and, more often than not, its results are very difficult to apply in 
practice or use as a guide. I will therefore settle here on a treatment that captures 
the essential physics of the problem, provides the formulas that describe the beha- 
vior of thermal conductivity in metallic systems, and has some predictive power 
regarding the choice of materials in various applications where heat conduction is 
an important issue. 

Because of the obvious practical relevance of noble metals such as copper, gold 
and silver and of various alloys containing copper, there has always been a strong 
interest in comparing the current and heat conducting characteristics of these ma- 
terials. It is therefore no surprise that certain important empirical relations were 
discovered nearly 50 years before the concept of an electron was firmly established 
and some 80 years before the notion of band structure emerged from the quantum 
theory of solids. Such is the case of the Wiedemann-Franz law, 6 which states that, 
at a given temperature, the thermal conductivity of a reasonably pure metal is 
directly related to its electrical conductivity; in other words, by making a simple 
measurement of electrical resistivity (conductivity a equals inverse resistivity f> '), 
we essentially know how good such a metal will be as a heat conductor, 

-= up = LT. (2) 

<7 

In many respects, the Wiedemann-Franz law and the constant L that relates the 
thermal conductivity of pure metals to their electrical conductivity, called the Lor- 
enz number, 7 are the fundamental tenets in the theory of heat conduction in metals. 
Most of our discussion will focus on the conditions under which this law is valid 
and on deviations that might arise. Such knowledge is of great interest: on one 
hand it allows a fairly accurate estimate of the electronic thermal conductivity 
without actually performing the rather tedious thermal conductivity measurements ; 
on the other hand, deviations from the Wiedemann-Franz law inform on the 
presence and strength of particular (inelastic) scattering processes that might influ- 
ence the carrier dynamics. 

I wish to stress right here that the Wiedemann-Franz law strictly compares the 
electronic thermal conductivity with the electrical conductivity. In metals that have a 
substantial phonon contribution to the overall thermal conductivity (impure metals 
and alloys), forming the ratio with the measured values of the thermal and electrical 
conductivities will naturally lead to a larger magnitude of the constant L because 
the overall thermal conductivity will contain a significant contribution from pho- 



nons. 
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1.2. The Drude Model 

The first important attempt in the theoretical understanding of transport processes 
in metals is due to Paul Drude. 8 This comes merely three years after J. J. Thomson 9 
discovered the electron, an elementary particle that clearly fits the role of a carrier 
of charge in metals regardless of the fact that the convention for the positive 
current direction would much prefer a positively charged carrier. 

The classical free-electron model of metals developed by Drude builds on the 
existence of electrons as freely moving noninteracting particles that navigate 
through a positively charged background formed by much heavier and immobile 
particles. At his time, Drude had no idea of the shell structure of atoms — nowadays 
we might be more specific and say that the gas of electrons is formed from the 
conduction or valence electrons that are stripped from atoms upon the formation of 
a solid while the positively charged ions (nuclei surrounded by the core electrons) 
are the immobile particles located at the lattice sites. In any case, the classical 
electron gas is then described using the language of kinetic theory. In this picture, 
the straight-line thermal motion of electrons is interrupted by collisions with the 
lattice ions. Collisions are instantaneous events in which electrons abruptly change 
their velocity and completely “forget” what their direction of motion was just prior 
to the collision. Moreover, it is assumed that thermalization occurs only in the 
process of collisions; i.e., following a collision, although having a completely ran- 
dom direction of its velocity, the speed of the electron corresponds to the tempera- 
ture of the region where the collision occurred. 

In defining conductivities, it is advantageous first to introduce the respective 
current densities. Since electrons carry both charge and energy, their flow implies 
both electric and heat currents. The term electric current density , designated J e , is 
understood to represent the mean electric charge crossing a unit area perpendicular 
to the direction of flow per unit time. Let us assume a current flows in a wire of 
cross-sectional area A as a result of voltage applied along the wire. If n is the 
number of electrons per unit volume and all move with the same average velocity 
v, then in time dt the charge crossing the area A is -nevAdt. Hence, the current 
density is 

J e = —nev. (3) 



In the absence of driving forces (electric field or thermal gradient), all directions of 
electron velocity are equally likely and the average velocity is zero; hence there is 
no flow of charge or energy. 

If an electric field e is applied to an electron gas, electrons accelerate between the 
collisions and acquire a mean (drift) velocity directed opposite to the field: 



dv es 

dt m 



( 4 ) 



Integrating Eq. (4), one obtains 

PF 

v(i) = --i + v( 0) 



( 5 ) 



where v(0) is the velocity at t-0, i.e., right after the collision at which time the 
electric field has not yet exerted its influence. The apparently unlimited rise in the 
velocity v(t) with time is checked by a subsequent collision in which the electron is 
restored to a local thermal equilibrium. There will be a range of times between 




Sec. 1 • INTRODUCTION 



25 



collisions, so to find out the mean value of the velocity over all possible times t 
between collisions, we need to know how these times are distributed. The prob- 
ability that an electron which has not collided with an ion during the time t will 
now suffer a collision in the time interval t + At is 

p-t/r 

P(t)dt = dt. (6) 

T 



The mean velocity then follows from 

ea ~ ea 

v= t = 

m in 



e f / T er 

t dt = e; 

r m 



(7) 



i.e., the mean time t between collisions is equal to r, the parameter frequently called 
the collision time or relaxation time. After Eq. (7) is substituted into Eq. (3), the 
electric current density becomes 



r nrar 

Je = = CT£, 

m 



( 8 ) 



from which one obtains the familiar expression for the electrical conductivity: 

ne 2 r 

a = . 

m 



(9) 



Thus, the Drude Theory predicts the correct functional form for Ohm’s law. By 
defining the mean free path between the collisions l e as 



l e = VT, 


(10) 


We can write Eq. (9) as 




ne 2 l e 

a = . 


(11) 


mv 


In Drude’s time, the mean speed of electrons was 
tion theorem, namely 


calculated by using the equiparti- 


^m(v) 2 = ^ k B T . 


(12) 



Thus, with the obtained mean velocity and using the experimental values of the 
electrical conductivity, the mean free path of metals at room temperature invariably 
came out in the range 1-5A, seemingly providing an excellent support for Drude's 
model that assumed frequent collisions of electrons with ions on the lattice sites. 

Similar to the electric current density, one can define the thermal current density 
Jq as a vector parallel to the direction of heat flow with a magnitude equal to the 
mean thermal energy per unit time that crosses a unit area perpendicular to the 
flow. Because the speed of an electron relates to the temperature of the place where 
the electron suffered the most recent collision, the hotter the place of collision the 
more energetic the electron. Thus, electrons arriving at a given point from the 
hotter region of the sample will have higher energy than those arriving at the 
same point from the lower-temperature region. Hence, under the influence of the 
thermal gradient there will develop a net flux of energy from the higher-temperature 
end to the lower-temperature side. If we take nil as both the number of electrons 
per unit volume arriving from the higher-temperature side and the same density of 
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electrons arriving from the lower-temperature region, and assuming that the ther- 
mal gradient is small (i.e., the temperature change over a distance equal to the 
collision length is negligible), it is easy to show that the kinetic theory yields for 
the heat current density the expression 

Jq = ^c„(-VT) = - Ke VT. (13) 



Here l e is the mean free path of electrons, c v is their electronic specific heat per unit 
volume, and tz, is the thermal conductivity of electrons. Equation (13) is a state- 
ment of the Fourier law with the electronic thermal conductivity given as 




(14) 



Relying again on the classical description of the electron gas, Drude writes the 
specific heat of electrons per unit volume c v in terms of the molar specific heat c m : 



N a 



3 

2 



NaIzb 





(15) 



Substituting into Eq. (14) yields the classical expression for the thermal conductiv- 
ity of electrons 



1 „ _ 3 nk B i e v 

K e = -zteV-ksn = — - — . 

3 2 2 



(16) 



By forming the ratio nja, one arrives at the Wiedemann-Franz law: 

Ke _ nk B £ e v/2 _ k B m(v) 2 _ 3 / 2 
a ne 2 £ e /mv 2e 2 2\eJ 



In the last step in Eq. (17) the mean square velocity is taken in its Maxwell- 
Boltzmann form 3 k B T/m. 

Equation (17) is a crowning achievement of the Drude theory — it accounts for an 
empirically-well-established relation between thermal and electrical conductivities of 
metals, and the numerical factor | ( k B /e ) 2 = 1.1 x 1(T 8 V 2 /K 2 is only a factor of 2 or 
so smaller than the experimental data . 

Overall, the success of the Drude theory appears impressive, especially when 
viewed from the perspective of the early years of the 1900s, and for the next three 
decades it formed the basis of understanding of the transport properties of metals. 
In retrospect, the success of this theory is quite fortuitous and stems from a lucky 
cancellation of a couple of large and grossly erroneous parameters that have their 
origin in the use of classical statistics. The most blatant error comes from the 
classically calculated value of the electronic specific heat, which grossly overesti- 
mates the actual electronic contribution to the heat capacity of metals. With its 
value of 3R/2 per mole, the classical electronic term is comparable to that due to 
the lattice and would result in far too large a specific heat of metals. Experimen- 
tally, no such classical electronic contribution was ever seen. This outstanding 
puzzle in the classical theory of metals lasted until it was recognized that electrons 



In his original paper, Drude makes a mistake in his evaluation of the electrical conductivity, which 
comes out only one-half of what Eq. (11) implies. With such an erroneous result his value for the Lorenz 
constant is 2.2 xl0~ 8 V 2 /K 2 in excellent agreement with experimental values. 
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are fermions and their properties must be accounted for with Fermi-Dirac statis- 
tics. 

In the expression for the Wiedemann-Franz law, Eq. (17), this large error is 
compensated by a comparably large error arising from the mean square electronic 
speed which, in the classical Maxwell-Boltzmann form, is a factor of 100 smaller 
than the actual (Fermi) velocity of electrons. The fortuitous cancellation of these 
two large errors leaves the Lorenz constant approximately correct except for its 
numerical factor of 3/2. With Fermi-Dirac statistics this numerical factor becomes 
equal to 7t 2 /3, in which case it is often stated that the Lorenz constant assumes its 
Sommerfeld value: 

To = -y (~) 2 — 2.44 x 1CT 8 V 2 K~ 2 . (18) 



Such fortuitous cancellation does not extend to the electrical conductivity in Eq. 
(11) and classical statistics is responsible for underestimating the mean free path by 
a factor of 100. In reality, l e should be on the order of several hundred angstroms 
rather than a few angstroms, as obtained by Drude. Actually there is another 
reason why Drude’s viewpoint of very frequent scattering of electrons on ions is 
not a good physical picture. With the advent of quantum mechanics it was soon 
recognized that, because the electrons have a wave nature (apart from their cor- 
puscular character), a perfectly periodic lattice presents no resistance to the current 
of electrons and the electrical conductivity in that case should be infinite. It is only 
because of deviations from perfect periodicity that the electrical conductivity attains 
its finite value. Imperfections can be of either dynamic nature and intrinsic to the 
structure, such as displacements of ions about the lattice sites due to thermal 
vibrations, or they can have a static character, such as impurities, vacancies, inter- 
stitials, and other structural defects. 

Another less than satisfactory outcome of the Drude theory that is a direct 
consequence of the classical treatment concerns the predicted temperature depen- 
dence of electrical resistivity. With the Maxwell-Boltzmann mean electron velocity, 
v = (3 ksT/m) 1 ^ 2 , the temperature dependence of resistivity is expected to be pro- 
portional to 



1 = jny ~ 1/2 

a ne 2 £ e 



(19) 



Experimentally, the resistivity of metals at temperatures above the liquid nitrogen 
temperature is typically a linear function of temperature. Again, the discrepancy 
with the classical Drude model is eliminated when quantum mechanics is used, 
namely when the mean velocity is taken as the Fermi velocity. 

The fact that we talk about the Drude model some 100 years after its inception 
indicates that, in spite of its shortcomings, the model provided convenient and 
compact expressions for electrical and thermal conductivities of metals and, with 
properly calculated parameters, offered a useful measure of comparison between 
important transport properties. Table 1 presents thermal and electrical conductiv- 
ities of metals together with their Lorenz ratio, all referring to a temperature of 273 
K. 
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TABLE 1 Thermal Conductivity of Pure Metals at 273 K.“ 



Metal 


k (W/m-K) 


p (10 -8 Sim) 


L (10~ 8 V 2 /K 2 ) 


Ref. 


Ag 


436 


1.47 


2.34 


11 


A1 


237 


2.43 


2.10 


12 


Au 


318 


2.03 


2.39 


13 


Ba 


23.3 


29.8 


2.55 


14 


Be 


230 


2.8 


2.36 


15 


Ca 


186 


3.08 


2.13 


16 


Cd 


100 


6.80 


2.49 


15 


Ce 


11.2 (291 K) 


80.0 


~ 3.5 


17 


Co 


99 (300 K) 


5.99 


1.98 (300 K) 


18 


Cr 


95.7 (280 K) 


11.8 


4.11 


19 


Cs 


37 (295 K) 


18.0 


2.5 (295 K) 


20 


Cu 


402 (300 K) 


1.73 (300 K) 


2.31 (300 K) 


21 


Dy 


10.4 (291 K) 


93 


3.75 (291 K) 


22 


Er 


13.8 (291 K) 


79 


3.75 (291 K) 


22 


Fe 


80.2 (280 K) 


8.64 


2.57 (280 K) 


23 


Ga (|| c) 


16.0 


50.3 


2.95 


19 


(II a) 


41.0 


16.1 


2.41 


19 


(II b) 


88.6 


7.5 


2.43 


24 


Gd 


9.1 (291 K) 


128 


4.2 (291 K) 


22 


Hf 


22.4 (293 K) 


31.0 (293 K) 


2.45 (293 K) 


25 


Hg (||) 


34.1 (197 K) 


14.6 (197 K) 


2.53 (197 K) 


26 


(A) 


25.9 (197 K) 


19.3 (197 K) 


2.55 (197 K) 


26 


Ho 


11.8 (300 K) 


78.0 (300 K) 


3.2 (300 K) 


27 


In 


81.0 (280 K) 


8.25 (280 K) 


2.39 (280 K) 


28 


Ir 


149 (277 K) 


4.70 


2.57 (277 K) 


29 


K 


98.5 


6.20 


2.24 


30 


La 


14.0 (291 K) 


59 


2.9 


22 


Li 


65 


8.5 


2.05 


31 


Lu 


16.2 (291 K) 


~ 50 


3.3 (291 K) 


22 


Mg 


153 (301 K) 


4.5 (301 K) 


2.29 (301 K) 


32 


Mn (a) 


7.8 (291 K) 


137 


4.0 (291 K) 


22 


Mo 


143 


4.88 


2.56 


33 


Na 


142 


4.29 


2.23 


34 


Nb 


51.8 (280 K) 


13.3 


2.53 (280 K) 


35 


Nd 


16.5 (291 K) 


58 


3.7 (291 K) 


22 


Ni 


93 (280 K) 


6.24 


2.19 (280 K) 


36 


Os 


87 (323 K) 


8.3 


2.7 (323 K) 


29 


Pb 


35.5 


19.2 


2.50 


37 


Pd 


71.7 


9.74 


2.57 


38 


Pr 


12.8 


65 


3.1 (280 K) 


39 


Pt 


71.9 (280 K) 


9.82 


2.59 (280 K) 


23 


Pu 


5.2 (298 K) 


~ 130 


2.48 (298 K) 


40 


Rb 


55.8 


11.3 


2.30 


41 


Re 


49 


16.9 


3.05 


42 


Rh 


153 (280 K) 


4.35 


2.46 (280 K) 


43 


Ru 


110 (280 K) 


6.7 


2.72 (280 K) 


43 


Sc 


21.8 


44 


4.3 


44 


Sm 


13.4 (291 K) 


90 


4.3 (291 K) 


22 


Sn 


64 


10.6 


2.48 


15 


Sr 


51.9 


11.0 


2.18 


14 


Ta 


57.7 (280 K) 


12.1 


2.56 (280 K) 


45 


Tb 


10.4 (291 K) 


110 


4.25 (291 K) 


22 


Tc 


51 (300 K) 


16.7 


- 3.4 (300 K) 


46 
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Metal 


k (W/m-K) 


p (10- 8 Sim) 


L (10“ 8 V 2 /K 2 ) 


Ref. 


Th 


49.3 


13.9 


2.56 


47 


Ti 


22.3 


40 


3.25 


48 


T1 


50.6 


15 


2.8 


49 


U 


28 (278 K) 


24 


2.8 (278 K) 


50 


V 


35 (260 K) 


18.9 


2.4 (260 K) 


51 


w 


183 (280 K) 


4.85 


3.27 (280 K) 


23 


Y 


15.9 (291 K) 


~ 52 


2.9 (291 K) 


22 


Zn 


114.5 (283 K) 


5.5 


2.31 (283 K) 


52 


Zr 


20.5 (323 K) 


39 


3.4 (323 K) 


53 


As (pc) 


27 (298 K) 


31.7 (298 K) 


2.87 (298 K) 


54 


(J_ c-axis) 


51.0 (300 K) 


28.4 (300 K) 


4.83 (300 K) 


55 


(|| c-axis) 


29.0 (300 K) 


38.6 (300 K) 


3.73 (300 K) 


55 


Bi (pc) 


7.8 (298 K) 


148 (298 K) 


3.87 (298 K) 


56 


(_L c-axis) 


9.8 (300 K) 


112 (300 K) 


3.65 (300 K) 


57 


(|| c-axis) 


6.0 (300 K) 


135 (300 K) 


2.69 (300 K) 


57 


Sb 


18.2 


43 


2.87 


58 



“ The data also include values of electrical resistivity and the Lorenz number. Most of the data are taken 
from the entries in the extensive tables in Landolt-Bornstein, New Series 10 . 



2. SPECIFIC HEAT OF METALS 



One of the major problems and an outstanding puzzle in the Drude picture of 
metals was a much too large contribution electrons were supposed to provide 
toward the specific heat of metals. While the classical treatment invariably sug- 
gested an additional 3i?/2 contribution to the molar specific heat of metals (on 
top of the 3 R term due to the lattice vibrations), the experimental evidence was 
unequivocal: the specific heat of simple solids at and above room temperature did 
not distinguish between insulators and metals, and the measurements yielded values 
compatible with the Dulong-Petit law; i.e., the molar specific heat of 3R = 6 cal 
mole -1 K -1 = 25 J mole 'K This puzzle was solved when electrons were treated 
as quantum mechanical objects and described in terms of Fermi-Dirac statistics. 
We briefly outline here the key steps. 

Having half-integral spin, electrons are classified as Fermi particles (fermions) 
and are subject to the Pauli exclusion principle. As intrinsically indistinguishable 
particles, symmetry requirements imposed on the wave function of electrons are 
such that no two electrons are allowed to share the same quantum state, i.e., to 
possess identical sets of quantum numbers. This fundamental restriction has far- 
reaching consequences as to how to treat the electrons, how to “build up” a metal, 
and, for that matter, how to assemble elements in the periodic table. The statistics 
that guarantee that the Pauli principle is obeyed were developed by Fermi and 
Dirac, and the distribution of fermions governed by these statistics is referred to 
as the Fermi-Dirac distribution f(E, T ) : 

/(-£> T ) = e ( E -o/k B T + i ' ( 20 ) 

Here ( is the chemical potential, often referred to as the Fermi level or Fermi 
energy and written as E F instead of (j. Equation (20) expresses the probability of 
finding an electron in a particular energy state E and has a very interesting prop- 
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erty. At T = 0 K, the function f(E,T ) is a step function that sharply divides the 
occupied states from the unoccupied states, with the Fermi energy playing the 
demarcation line between the two. At E < E F , all available single-electron states 
are fully occupied while for E > E F all states are empty. The Fermi energy is thus 
the highest occupied energy state. As the temperature increases, a certain degree of 
“smearing” takes place at the interface between the occupied and unoccupied 
states, with states close to but below E F now having a finite probability of being 
unoccupied, while states close to but slightly above E F have a finite chance to be 
occupied. The range of energies over which the distribution is smeared out is fairly 
narrow and amounts to only about 4 k F T around the Fermi level. Since the Fermi 
energy of typical metals is several electron-volts (eV) while the value of 4 k F T at 300 
K is only about 100 meV, the smearing of the distribution is indeed small on the 
fundamental energy scale of metals. Nevertheless, this lack of sharpness in the 
distribution of electrons is all-important for their physical properties. Because the 
deep lying electron states are fully occupied with no empty states in the vicinity, 
electrons occupying such states cannot gain energy nor can they respond to external 
stimuli such as an electric field or thermal gradient. Only electrons within the 
smeared-out region of the distribution have empty states accessible to them and 
can thus absorb the energy. Of the entire number of free electrons constituting a 
metal, it is thus only a fraction 4 ksT/Ep (~0.01 ) of them that are involved in the 
specific heat. This is the reason why one does not detect electronic contribution to 
the specific heat of metals at and above ambient temperatures; the contribution of 
electrons is simply too small on the scale of the Dulong-Petit law. 

To calculate the specific heat due to electrons exactly within the single-electron 
approximation, apart from their distribution function, one also needs to know the 
number of energy states per unit volume in a unit energy interval, i.e., the density of 
states, D{E). This quantity is model dependent. Regardless of the model chosen for 
D{E), its product with the distribution function integrated over all energies of the 
system must yield the density of electrons per unit volume: 

OO 

n = J D(E)f(E, T)dE. (21) 

0 

When electrons are heated from absolute zero temperature to some finite tempera- 
ture T . , they absorb energy and increase their internal energy by 

oo 

U(T) = J ED(E)f{E,T)dE - J ED(E)dE. (22) 

o o 

The second integral in Eq. (22) represents the internal energy of electrons U{ 0) at 
7=0 K. Electrons have this non-zero internal energy at absolute zero purely be- 
cause they obey the Pauli exclusion principle, and only two electrons out of n 
electrons (one with spin up and one with spin down) can occupy the lowest energy 
state. All other electrons have to be accommodated on the progressively higher 
energy levels, hence a significant internal energy even at T - 0 K. By differentiating 
with respect to temperature and employing a trick of subtracting a term the value 
of which is zero but the presence of which allows one to express the result in a 
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convenient form,* we obtain 

OO 

Cv = J (E-E f )D(E) df ^ T T) dE. (23) 

o 

Assuming that the density of states does not vary much over the narrow tempera- 
ture range where the distribution function is smeared out around the Fermi energy, 
(i.e., taking the density of states at its value of the Fermi energy E F and carrying 
out the partial derivative of the distribution function), one obtains 

+oo 

/ 2 x 2 

, X ~ e \ 2 dx = \ klTD(E F ) . (24) 

(e x + 1) j 

— OO 

This is a general result for the specific heat per unit volume of electrons, indepen- 
dent of a particular form of the density of states. For a free-electron gas, it is easy 
to show that the density of states becomes 

D(E) = [ E l/2 . (25) 

V ' 2 tt 2 ^ 3 V ' 

In that case, the density of electrons per unit volume can be conveniently expressed 
as 

n = ^D(E F )E F , (26) 

and with the definition of the Fermi temperature Tp = E F /k B , the heat capacity per 
unit volume becomes 

Cv=\nk B ^r- (27) 

Comparing with Eq. (15), one sees that the effect of Fermi-Dirac statistics is to 
lower the specific heat of electrons by a factor of (tt 2 /3)(T /TV). Because the Fermi 
temperature of metals is very high, ~10 4 -10 5 K, this factor is less than 0.01 even at 
room temperature. 

To get from the specific heat per unit volume into a more practically useful 
quantity — specific heat per mole of the substance — Eq. (27) must be multiplied 
by the volume per mole, ZNa/ti, where ZNa is the number of electrons in a 
mole of metal of valence Z and N A stands for Avogadro’s number: 

C = ~2 Zk B N A ^p = = 7T. (28) 

In Eq. (28), R is the universal gas constant and 7, the coefficient of the linear 
specific heat of metals, is often called the Sommerfeld constant. Although the 
magnitude of the electronic specific heat is small, and certainly at ambient tempera- 
ture it is negligible compared to the lattice contribution, it gains importance at low 



Because of the high degeneracy of electrons, the electron density n in Eq. (21) is a constant and by 
multiplying it by the value of the Fermi energy, the product EFn remains temperature independent. 
Differentiating it with respect to temperature gives 0 = J EpD(E){FIT)dE. This term is then subtracted 
from the derivative of Eq. (22) with respect to temperature to obtain Eq. (23). 
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temperatures. Because the lattice specific heat decreases with a much faster T 3 power 
law, there will be a crossover temperature (typically a few kelvins) below which a 
slower, linear dependence of the electronic contribution will start to be manifested 
and eventually will become the dominant contribution. 

It is interesting to compare experimental values of the coefficient 7 of selected 
metals with the corresponding free-electron values given by Eq. (28). In essence, 
such a comparison (see Table 2) informs how appropriate it is to describe the 
electrons of each respective metal as being free electrons. We note that for the 
noble metals and the alkaline metals the free-electron description is a reasonable 
starting point, while for most of the transition metals the free-electron picture is 
grossly inadequate. 

The large deviations between the experimental and free-electron 7-values ob- 
served in the case of transition metals are believed to arise from a very high density 
of states of the (/-electrons near the Fermi energy. Transition metals are character- 
ized by their incomplete (/-shells. When transition metals form by bringing together 
atoms, the two outer s-electrons split off and go into a wide conduction band. The 
more tightly bound (/-electrons also form a band, but this is a rather narrow band 
because the overlap between the neighboring (/-orbitals is small. Nevertheless, this 
band must be able to accommodate 10 electrons per atom, hence its large density of 
states. On the other hand, semimetals such as Sb and Bi possess very low values of 
the 7 parameter because their carrier density is several orders of magnitude less 
than that of a typical metal. 

3. THE BOLTZMANN EQUATION 

When one talks about transport phenomena, one means processes such as the flow 



TABLE 2 Comparison of Experimental and Calculated (free-electron) Values of the 7 Coefficient of the 
Electronic Specific Heat for Selected Metals."' 



Metal 


Valence 


Free-electron 7 
(10~ 3 J mokr'K- 2 ) 


Experimental 7 
(10 3 J mole-' K~ 2 ) 


Na 


1 


1.09 


1.7 


K 


1 


1.67 


1.7 


Ag 


1 


0.64 


0.66 


Au 


1 


0.64 


0.67 


Cu 


1 


0.50 


0.69 


Ba 


2 


1.95 


2.72 


Ca 


2 


1.51 


2.72 


Sr 


2 


1.80 


3.64 


Co 


2 


0.61 


4.98 


Fe 


2 


0.64 


5.02 


Ni 


2 


0.61 


7.02 


A1 


3 


0.91 


1.30 


Ga 


3 


1.02 


0.63 


Sn 


4 


1.39 


1.84 


Pb 


4 


1.50 


2.93 


As 


5 


1.30 


- 


Sb 


5 


1.63 


0.628 


Bi 


5 


1.79 


0.08 



a. Experimental values are taken from Ref. 2. 
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of charge or flow of heat in solids. One is usually concerned about the steady-state 
flow, i.e., the flow established as a result of the influence of external driving forces 
(electric field or thermal gradient in our case) and the internal scattering processes 
tending to restore the system to equilibrium. A steady state therefore must be 
distinguished from the equilibrium state, and transport theory is thus a special 
branch of nonequilibrium statistical mechanics. The deviations from the truly equi- 
librium state, however, are usually small. The questions are: how are the driving 
forces and the scattering processes interrelated? and how does the population of the 
species under investigation evolve as a function of time? The answers are provided 
by the Boltzmann equation. Although there are other approaches describing how to 
treat the problem and arrive at formulas for the transport parameters, the Boltz- 
mann equation has proven itself time and again as the most versatile approach and 
the one that is reasonably easy to track. Moreover, it yields a form of the expres- 
sions for the transport parameters that are intuitive and that can readily be com- 
pared with the experimental results. In this section we show how one uses the 
Boltzmann equation to arrive at useful formulas for the transport parameters of 
metals. 

Metals are crystalline solids. The periodic potential of ions located at the crystal 
lattice sites is an essential feature of the system and must be incorporated in a 
realistic description of metals. Lattice periodicity has the following consequences: 
the momentum of a free electron is replaced by a quantity //k, where k is the Bloch 
wave vector of an electron. The wave vector contains information about the trans- 
lational symmetry of the lattice. The electron velocity, which in the free-electron 
model is simply v = plm , is replaced by the group velocity (l//;)<9E(k)/<9k, where E(k) 
is the Bloch energy of an electron. Wave functions, which in the free-electron 
picture were simply plane waves exp(/ kr ), must now reflect the periodicity of the 
lattice. This is accomplished by introducing Bloch waves ip k ( r) = «t(r)exp(/ k . r ), i.e., 
plane waves modulated by a function u k ( r) that has the periodicity of the lattice. 

The fundamental issue concerns the distribution function of electrons in a metal, 

i.e., how the occupation number of electrons changes as a result of the influence of 
an electric field and thermal gradient, and the effect of various scattering processes 
that the electrons undergo. The equilibrium distribution of electrons is governed by 
the Fermi-Dirac function [Eq. (20)], which is independent of the spatial coordinate 
r because of the assumed homogeneity (in equilibrium, the temperature is the same 
everywhere). Away from equilibrium, the distribution may depend on the spatial 
coordinate r, since we assume that local equilibrium extends only over the region 
larger than the atomic dimensions, and of course on the time t. Moreover, because 
the energy of electrons is a function of the wave vector, we also assume that the 
distribution function depends on k. We therefore consider explicitly /(r,k,t). 
Temporal changes in /(r.kj) arise because of three influences: 

1 . Electrons may move in or out of the region in the vicinity of point r as a result of 
diffusion. 

2. Electrons are acted upon by driving forces such as an electric field or thermal 
gradient. 

3. Electrons are deflected into and out of the region near point r because they 
scatter on lattice vibrations or on crystal imperfections. This is the scattering 
influence on the distribution function, and the respective partial derivative is 
called the collision term. 
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The first two influences on /(r,k,t) have their origin in the Liouville theorem 
concerning the invariance of the volume occupied in phase space. Thus, the number 
of electrons in the neighborhood of point r at time t must equal the number of 
electrons in the neighborhood of point r - v(k)d? at the earlier time t - dt. Since the 
electrons are also subjected to an external electric field E, they accelerate and in 
time dt change their momentum /;k by (dkldt)dt = e zdt/h. (We assume that there is 
no magnetic field present.) Therefore, in analogy with the Liouville theorem for the 
spatial coordinate r, we consider the same volume invariance but this time for 
L-space; i.e., the electrons that are at point k at time t must have been at a location 
k - (~e) zdt/h at time t - dt. In the absence of any kind of collisions, we therefore 
can write 



/(r,M) = / 



. , , dt 

r — v(k)af, k+ ee — ,t — dt 



(29) 



However, collisions cannot be neglected since they change the population of elec- 
trons. Some electrons will be deflected away from the stream of electrons during the 
time interval dt as the electrons proceed from (r - \(t)dt, k + ezdt/h) to (r, k). 
However, it is also possible that some electrons that arrive at (r, k) at time t do not 
belong to the stream of electrons we considered at (r -\{t)dt, k + ezdt/h ) at time 
t - dt because they were deflected into the stream by collisions in the neighboring 
regions during the time dt. It is customary to capture both the “out-deflected” and 
“in-deflected” electrons in a term ( dfldt) coa . Thus, in time dt there will be a change 
of population due to scattering of magnitude (dfldt) co \\dt. Putting it all together, we 
write 



/( r, k ,t) = f 



. , , dt 

r — v(k)af , k + ez — , t — dt 
n 




(30) 



Expanding the equation to terms linear in dt, one obtains the Boltzmann equation : 



Of 

dt 



+ v • V r / 




(31) 



Equation (31) describes a steady state, not the equilibrium state. The terms on the 
left-hand side are frequently called the streaming terms. In principle, the collision 
term on the right-hand side contains all the information about the nature of the 
scattering. Using the quantum mechanical probability for transitions between the 
Bloch states ipk and i/V, specifically Ww oc | <k'|f?'|k> | 2 , where H' is the pertur- 
bation Hamiltonian, one can write the collision term as 



f9/(k)] 


1 _ v r 


{ dt ) 


' coll (27 t) 3 J 



dk'{[l - f(k)]W w f(k') - [1 - /(k')]WW(k)} 



(32) 



The two terms in braces represent transitions from the occupied state k' to an 
empty state [1 - /(k)J and a corresponding transition from the occupied state k 
to an empty state [1 - /(k')]. These terms are present to comply with the Pauli 
exclusion principle. Substituting Eq. (32) into Eq. (31), one gets a complicated 
integrodifferential equation for the distribution function /(r.kj). 

One of the most successful approaches for solving the Boltzmann equation relies 
on the use of a relaxation time approximation. The essence of this approach is the 
assumption that scattering processes can be described by a relaxation time r(k) that 
specifies how the system returns to equilibrium, i.e., how the distribution function 
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/(r,k,f) approaches its equilibrium value /°(k). One therefore writes the collision 
term as 



(df) /(k)-/°(k) _ ff(k) 
l dt ) con r(k) r(k) ’ 



where g(k) is the deviation of the distribution function /(k) from its equilibrium 
value /°(k), and r(k) is the relaxation time. Further considerable simplification of 
the transport problem is achieved by linearization of the Boltzmann equation. This 
is done simply by replacing the steady-state electron distribution in the gradients 
V r /(r,k,t) and V/,/(r.kj) by the equilibrium distribution function, i.e., by taking 
V r /(r,k,r) = V r /°(r,k) and VA./(r,k,?) = VA-/°(r,k). This seems eminently reasonable 
provided that the external fields are not too strong so that the steady-state distribu- 
tion is not too far from equilibrium. We can easily evaluate these two respective 
gradients (including specifically the dependence of the Fermi energy on tempera- 
ture) and rewrite Eq. (31) as 




E — Ep 
T 



V r T 




(34) 



This is the linearized Boltzmann equation. After some algebra, we can write, for the 
perturbed distribution function f(r,k,t), 

/(k) = /°(k)+ ( — ^)v(k)T(k)|e[e-^V r B F ] + (— V r T) j> . (35) 



This is a general form of the distribution function describing electron population 
perturbed by a weak electric field and a small temperature gradient. Before we 
proceed further it is important to comment on the quantity we call the relaxation 
time r(k). 

It is a major consideration whether r(k) is a true relaxation time in the sense of 
representing something like the time between collisions and, therefore, whether it 
faithfully represents the way the equilibrium state is approached. We can, of course, 
always formally introduce the parameter r(k) for any process, but if this quantity 
turns out to depend on a particular form of the field that created the out-of- 
equilibrium state, then we would need different functional forms of r(k) for differ- 
ent external fields and the meaning we associate with the relaxation time would be 
lost. Thus, whether the true relaxation time exists (i.e., whether the relaxation time 
approximation is valid) is a critical issue, and we shall see that it has a decisive 
influence on the discourse concerning the thermal conductivity at low temperatures. 



3.1. Transport Coefficients 

Knowing how a weak electric field and a thermal gradient perturb the population 
of electrons, we may now proceed to inquire what currents arise when the system is 
allowed to reach steady state. To do that, it is customary to define current densities. 
If we consider a flow of particles, we know that a volume element dk around point 
k in reciprocal space will contribute a flux of particles of magnitude ( 1/87T 3 ) y(k)c/k, 
where v(k) is the particle velocity. Since we consider electrons, there are two spin 
states for every Ar-state and the electrons are distributed according to the Fermi- 
Dirac function. Summing over all available states and assigning a charge e to each 
electron, we obtain the current density: 
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J e = J v(k)/(k)dk. 

Similarly for the heat current density, we write 

Jq = ^ 3 / v(k)[£7-C]/(k)dk. 



(36) 



(37) 



Equation (37) contains a term E — where £ is the chemical potential (free energy, 
which in metals can be taken as the Fermi energy) rather than just energy E. This is 
necessary because thermodynamics tells us that heat is internal energy minus free 
energy. 

Since equilibrium distributions do not generate currents (there is no spontaneous 
current flow), substituting the perturbed distribution function in Eq. (35) into Eq. 
(36) causes the integral containing the equilibrium distribution to vanish. The 
remaining term is 



Jr = 



2e 

8-7T 3 



e 2 r 

4ir 3 h 



er 

4x^h 



J v(k)v(k)r(-|Q [e£-V£ + ^(-VT) 
v(k)v(k) (~§s] [ £ ~j v< 

< k H k )(-D r£ - c ' 



dk 



n dE 



(- VT ) t\ dE ■ 



(38) 



Similarly, substituting Eq. (35) into Eq. (37) results in 

J « = 8? / f ) [« - VC + ^ (- VT) 



(A - C)dk 



= ^ / v ( k M k ) 



E--VC 

e 






T 

47 r 3 






(39) 



where, with the aid of dk = dSdk ± = dSdE /\X7 k E\, we converted an integral over a 
volume in A'-space into an integral over surfaces of constant energy. Defining now 
the integral 

* 7/voov*,,*- cr(-£)j>. 



Kn = 



4ir 3 h 



we can express Eqs. (38) and (39) in terms of K n : 

J e = e 2 A' 0 £ + ^(-VT) 



(40) 



(41) 



J Q = eAh£ + ^(-VT). 



(42) 
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Equations (41) and (42) reveal an important point — electric current J e can be 
generated not just by an external electric field but also by a thermal gradient. 
Likewise, heat current J q arises as a result of an electric field and a thermal 
gradient. In essence, there are interactions between the electric and heat currents. 

Integrals K n ; defined in Eq. (40), can be evaluated, using a general theorem for 
integrals of the Fermi function over energy that are of the form 

+oo +oo 

J *(. E)f(E)dE = j' $(J5) [~%)dE, (43) 



where ip(E) = dA>( E)ldE. Because -df°ldE has an appreciable value only within a 
few k B T of C, a smoothly varying function of energy $(£) can be expanded about 
E = (, where hopefully only the first few terms suffice for accuracy : 



$( 25 ) = <lKC) + £ 

n= 1 



\(E-C) n ] 


'd n <S>(E)' 


n\ 


dE n 



(44) 



The integral in Eq. (43) then becomes 




dE = 4>(C) 



7T 2 ( k B T ) 2 

6 



d 2 <S>{E) 



dE 2 



E= C 



(45) 



and only terms even in n contribute because df°ldE is an even function of E — (. 

Equations (41) and (42) serve to define transport coefficients in zero magnetic 
field. The usual isothermal electrical conductivity follows from a condition that the 
thermal gradient is zero: 

J e = e~K 0 z. (46) 



With the explicit form of K 0 , this can be written as 



Jr = 



e 2 t{E f ) 
47 r 3 h 



e=e f 



v(k)v(k) Fi)T E = <rs ’ 



(47) 



where a is the electrical conductivity. In a general crystal structure the current 
density need not be parallel to the electric field, and the electrical conductivity is 
a tensor, 



_ e 2 r(E F ) f' Vi(k.)vj(k)dS E 

' ij ~ 4 n 2 h J |v(k)| 

e=e f 



(48) 



For cubic metals and for isotropic (i.e., polycrystalline metals), the tensor reduces 
to a scalar. Placing now the electric field e along, say, the x-direction, we obtain, 
making use of v 2 = (l/3)u 2 : 



_ e 2 r(E F ) I" 

127 r 3 h J 

e=e f 



vdS F = 



127r 3 h 



irdSF 



e=e f 



(49) 



In this equation, we introduced the mean free path of electrons as | e = v.t (E f ). 

Equation (49) is the basic formula for the electrical conductivity. It states clearly 
that only electrons near the Fermi level can contribute to the transport process 
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since the integral is over the Fermi surface T(k) = Ep. This, of course, is in full 
agreement with the Pauli exclusion principle. 

One might think that the thermal conductivity is obtained simply by taking it as 
a coefficient of the thermal gradient in Eq. (42), i.e., under the conditions that an 
external electric field is zero. Elowever, this is not how thermal conductivity is 
measured. Although it is difficult to control electric fields, it is a straightforward 
matter to ensure that no electric current passes through the sample while the 
thermal conductivity is measured. We simply leave the sample in open-circuit con- 
figuration. Therefore, to obtain the desired coefficient of thermal conductivity, we 
set the current J e = 0 in Eq. (41) and express an electric field in terms of a thermal 
gradient : 



£ = 



eT 



VT. 



(50) 



Substituting Eq. (50) into Eq. (42) yields 

1 



Jq = J, [Ki - K\K~ x K\\ (— VT) = k(-VT). 



(51) 



The coefficient n is the electronic thermal conductivity. The term K\K~ l K\ in Eq. 
(51) is very small for metals and, for present purposes, can be neglected. The 
thermal conductivity is therefore 



n = ^[K 2 -K 1 K; l K l ] 



(52) 



If we define a function <&(£) = (E — E F ) 2 K 0 (E), it can be shown, using Eq. (45) and 
noting that the term [ ( I>(E)] B=Bf . = 0, that K 2 is related to K a via 



K 2 = AhA 2 K„(Er) = K 0 (E F ) 



(53) 



Substituting Eq. (53) into Eq. (52), we obtain the thermal conductivity as 

JLl2rr - 1 (k D \ 2 

Ta 



K 2 n~k~ B T 
K = ^r = — 3 



K 0 (E F ) = 



(?) 



(54) 



We recognize the ratio k/ct as the Eorenz number of the Wiedemann-Franz law. 
Equation (54) is to be compared to the Drude result [Eq. (18)], where the only 
difference is in the factor 7r 2 /3 = 3.29 in place of 3/2, arising because Eq. (54) is 
derived with the proper statistical description of the electrons. Eq. (54) is quite a 
general result except for the fact that scattering processes must be elastic. 

Now that the coefficients governing the flow of charge, a, and the flow of heat, k, 
have been established, two other coefficients remain to be determined based on Eqs. 
(41) and (42). These represent the interference terms between the electric and ther- 
mal fields. 

Let us assume that we set up a temperature gradient across a sample that is in 
open-circuit configuration; i.e., there is no electric current J e . This gradient gives 
rise to an electric field e: 



£ = 



K~ l K , 



VT = SVT. 



eT 



(55) 
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The coefficient S is the Seebeck coefficient, often, but somewhat unfortunately, 
called thermoelectric power: 



_ K'Ki 

eT 



(56) 



Using a relation between the integrals K\ and K a , 

dK 0 {E) 



Ah =y(fc B T ) 2 



dE 



J e=e f 



(57) 



one can write the Seebeck coefficient in the form 



7T 2 ek^T 


dK 0 (E)~ 




ikk^T 1 


'da(E)' 


7T 2 (k 


knT 


dlna(E) 


3 a 


dE 


E=E F 


3e (J 


dE 


E=E f 3 l ( 


1 


dE 



(58) 



The Seebeck coefficient is an important transport parameter and has a considerable 
bearing on technological applications such as thermocouple sensors and thermo- 
electric power conversion. The study of thermoelectric phenomena is an interesting 
topic but clearly beyond the purview of this chapter. A reader interested in thermo- 
electric effects might find it useful to consult monographs dedicated to this subject, 
among them Goldsmid, 59 Blatt et al., 60 and Tritt. 61 

The remaining term in Eqs. (41) and (42) is obtained by setting VT = 0, i.e., 
imposing isothermal conditions and relating the heat current to the electric current : 

J Q = Je = nj e , (59) 

where fl is called the Peltier coefficient. It relates heat generation or absorption in 
junctions of dissimilar metals to the electric current passed through the circuit. The 
Peltier effect should not be confused with Joule heating, which is a quadratic 
function of current and is always dissipative. The Peltier effect is linear in current 
and is at the heart of an important but niche technology of thermoelectric cooling. 
The materials of choice are, however, small-gap semiconductors rather than metals. 

For completeness, we also mention that it follows from Eqs. (56) and (59) that 
the Seebeck and Peltier effects are closely related, namely 

n = ST. (60) 



Having defined the transport coefficient and having established its form in terms 
of the transport integrals, we now address the key issue: the actual mechanism of 
transport. We will inquire about the dynamics of electrons and about the processes 
that limit both the electric and heat currents. Since the Wiedemann-Franz law plays 
such a pivotal role in the theory of transport and in the practical assessment of heat 
conducting ability of metals, we also inquire about the validity of this law in 
different regimes of transport. To address these issues, we rely on a quantum- 
mechanical description of the electron motion in a crystalline lattice and on the 
interactions of electrons with both intrinsic and extrinsic lattice imperfections. We 
also consider the influence of electrons interacting with each other. We should state 
right at the outset that a theoretical description of the transport mechanism is a 
difficult and challenging topic. As much as one would like to capture all nuances of 
the band structure of various metals and make the computations and analysis truly 
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metal specific, one soon finds out that the sheer complexity of the problem neces- 
sitates various approximations and simplifications that, in spite of valiant efforts, 
yield not much more than a qualitative solution. It is, indeed, unrealistic to expect 
that the theory will provide perfect agreement with the experimental data. What we 
hope for are correct predictions regarding the temperature dependence of various 
transport phenomena, capturing the trends among various classes of metals, and 
perhaps a factor of 2 consent regarding specific numerical values of the transport 
coefficients. Because electrical conductivity of metals is so closely tied to their 
electronic thermal conductivity, we consider its behavior first. This will serve as 
an excellent preparation for discussions of heat transport in metals. 



3.2. Electrical Conductivity 



In Eq. (49) we have a basic formula for electrical conductivity. To solve the integral 
in this equation, we must consider specific models. The most amenable one assumes 
that electrons are in a band that is strictly parabolic (nearly free electrons). In this 
case the velocity of electrons at the Fermi surface is simply v(E F ) = hk F /m * (m* 
being the effective mass), and the integral over the Fermi surface is 4nk 2 F . Collecting 
these results, we obtain the electrical conductivity 



e 2 r(E F ) 1 / hk F 
4n 3 h 3\ to* 




e 2 r(E F ) 0 _ e 2 r(E F ) ne 2 r(E F ) 

Q ? * 1 ? * ^ * 

37r z ra* 37 r z ra* rrv 



(61) 



In Eq. (61) we have recovered a result formally equivalent to the Drude formula, 
Eq. (9). This time, however, we are on a firm footing in terms of the theory. The 
relaxation time, which in Eq. (9) was purely a statistical quantity, is now a relaxa- 
tion time of the electrons on the Fermi surface. Only electrons that are within a 
thermal layer on the order of k F T of the Fermi energy can respond to an external 
electric field because they are the ones that can find empty states in their vicinity. 
Electrons deep down in the distribution have no empty states as neighbors and thus 
are unaffected. Metals have high conductivity because a relatively small number of 
them — those in the neighborhood of the Fermi energy — have very high velocities 
and not because all electrons drift “sluggishly” in response to the electric field. This 
is also easily understood from the following illustration that depicts how the out-of- 
equilibrium electron distribution actually looks. From Eq. (35), taking both gradi- 
ents V r T and V r E F equal to zero, we write 



m 



n k) 



er(k)ar^(k) ( ert 

n dE d k ' ''I h ) 



(62) 



where we used Taylor’s theorem to arrive at the second equality. Eq. (62) states that 
the steady-state distribution /( k) in reciprocal space (/c-space) is the same as the 
equilibrium distribution that is shifted by -(er//i)£. Therefore, depicting the equili- 
brium distribution of electrons by a solid curve in Fig. 1(a), we see that the effect of 
a constant electric field applied along the negative x-axis is to shift the entire 
distribution in the positive x-direction by ( er/h ) £. The displaced electron distribu- 
tion is indicated in Fig. la by a dashed curve. It is clear from the figure that only 
electrons near the Fermi edge are affected by the electric field ; electrons lying deep 
inside the distribution are completely ignorant of the presence of the field. If the 
external field is switched off, the steady-state distribution (the dashed curve) will try 
to relax back to the equilibrium form. The question is, how does the distribution in 
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e vt 

(a) (b) 



FIGURE 1 Schematic representation of the undisturbed (solid curves) and disturbed (dashed curves) 
Fermi distributions produced by (a) an electric field and (b) a temperature gradient. The overpopulated 
and underpopulated energy levels are marked with + and - signs, respectively. Solid and open circles in 
the upper panels represent the excess and deficiency of electrons relative to equilibrium distribution. An 
electric field shifts the entire distribution, whereas a temperature gradient creates asymmetry in the 
distribution function. The reader should note the distinction between the large-angle scattering (the 
so-called horizontal processes) and small-angle scattering but with a change in the electron's energy 
(the vertical processes). 



Fig. la relax? Obviously electrons must be taken from the region on the right and 
moved to the opposite side of the Fermi surface. Whether this can be achieved in a 
single interaction involving a large angular change or whether many smaller steps 
are needed to bring electrons back determines how effective is the scattering pro- 
cess. The average time needed to reestablish equilibrium is a measure of relaxation 
time r. This may turn out to be a different time than the average time between 
successive collisions, r 0 , and we must carefully distinguish these two times. Clearly, 
if the distribution can be relaxed in a single step, then r = r 0 . Flowever, if several 
collisions are needed to relax the distribution, the time r is larger than r 0 and is 
given by 

- = — (1 — (cos#)), (63) 

T T o 

where (cos <)) stands for an average over the scattering angle. The point is that from 
the perspective of electrical resistivity* it is the change in the electron velocity in the 
direction of the electric field which gives rise to resistivity. If each collision changes 
the direction by only a small amount, many collisions will be needed to accomplish 
the task of relaxing the distribution. Because the factor (1 - (cos 9 )) is sensitive to 
temperature, it will figure prominently in the discussion of electrical and thermal 
transport processes considered as a function of temperature. 

In a strongly degenerate system such as metals the carrier density is essentially 



Thermally perturbed distributions, in addition to a possibility of being relaxed via the mechanism 
shown in Fig. la, i.e., using the so-called horizontal processes, have an alternative and often more 
effective means of relaxation through very small angle but with a small change in energy, the so-called 
vertical processes illustrated in Fig. lb and discussed later. 
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temperature independent. It then follows from Eq. (61) that variations in the elec- 
trical resistivity with temperature come from the temperature dependence of the 
relaxation time r. With the aid of quantum mechanics, it is, in principle, possible to 
calculate all scattering probabilities and, therefore, to determine the relaxation time. 
However, Fermi surfaces of real metals are quite complicated. One either does not 
have available all the input parameters, or, if one has them, the theory is not 
capable of incorporating all the line nuances of the problem. In short, calculations 
of relaxation times are not a trivial task once one ventures beyond the comfort zone 
of spherical Fermi surfaces. 

We have already mentioned that it is somewhat arbitrary to write the collision 
term in the form given in Eq. (33), and that only in situations where r(k) is 
independent of the external fields is the relaxation time meaningful. Since in this 
chapter we only consider transport processes in zero magnetic field, we aim for a 
relaxation time that is common to both an electric field and a thermal gradient. 
Clearly, we would like to know when this is so; i.e., under what conditions does the 
relaxation time approximation yield a result equivalent to a rigorously calculated 
collision term in the Boltzmann equation? This is an important question that 
touches upon the validity of the Wiedemann-Franz law. It turns out that an 
electron population disturbed by an electric field and a thermal gradient will 
have the same relaxation time when electrons scatter elastically. Actually, as shown 
in Ref. 2, if the relaxation time approximation returns a solution for the distribu- 
tion function of the form 

/(k) = /°(k) + A(£)-k, (64) 

where a vector quantity A (E) depends on k only through its magnitude, then this 
solution is rigorous in the sense that the same result is obtained by solving the full 
Boltzmann equation. 

We must not forget that by using the relaxation time approximation we have 
merely introduced the relaxation time as a parameter but we have not solved the 
scattering problem; this can be done only by having a general solution to the 
Boltzmann equation. Tackling such a formidable challenge, one must resort to 
more sophisticated techniques. Of these, the variational principle has proved to 
be the most successful. Developed by Kohler, 62 it is based on the linearized Boltz- 
mann equation for an electric field and thermal gradient. It defines a function $(k) 
such that 

(65) 

The function $ is then expressed as a linear combination of appropriately chosen 
functions with variable coefficients. The aim is to optimize the coefficients so as to 
produce a trial function that most closely approximates the true solution. The 
Boltzmann equation is thus reduced to an extremum problem. The better the trial 
function, the better is the approximation. An interesting upshot of this treatment is 
its thermodynamic interpretation. One can view the effect of an external electric 
field as aligning the velocities of electrons and thus decreasing the entropy. On the 
other hand, scattering tends to disperse electrons, i.e., decrease the degree of order 
in the system and thus increase the entropy. The Boltzmann equation therefore 
expresses a relation between the decrease in entropy (by external fields) and the 
increase in entropy (by collisions). Since variational calculations test for the max- 
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imum, the solution of the Boltzmann equation coincides with the maximum entropy 
production in the collision processes. For more information on the variational 
principle as applied to the Boltzmann equation readers might wish to consult 
Refs. 3 and 63. 

For completeness, we mention that the Boltzmann equation is not the only 
plausible approach for the formulation of the transport problem. One may, for 
instance, use as a starting point linear response theory and derive the conductivity 
in terms of time correlation functions of current. This is the approach pioneered by 
Kubo 64 and by Greenwood: 65 

OO 

a = ( 66 ) 
0 



Other correlation functions (e.g., force-force correlations 66 ) have been tried to 
arrive at an expression for the conductivity. A different approach, based on the 
density matrix formalism, was developed by Kohn and Luttinger. 67 Even the op- 
tical theorem was explored to derive an expression for conductivity. 68 Although 
these approaches capture the a priori microscopic quantum nature of the problem, 
none of these delivers a clearly superior solution. Yes, there is considerable comfort 
and satisfaction that one can start from the fundamentally atomistic perspective 
and use the rigorous field-theoretical techniques to arrive at a solution. Flowever, 
the road to this solution is rather arduous, the treatment is less intuitive, and the 
results are often difficult to interpret. In spite of its shortcomings and a somewhat 
less than rigorous nature, the Boltzmann transport equation is not a bad starting 
point for the theoretical description of transport phenomena. 

From Eqs. (33) and (61) it follows that the central problem in transport theory is 
a calculation of the relaxation time r. In metals, three main scattering processes 
affect the electrical and thermal resistivities, (a) electrons can scatter on lattice 
defects such as foreign atoms (impurities) occupying the lattice sites; (b) electrons 
can be deflected via their interaction with lattice vibrations (phonons); and (c) 
electrons may interact with other electrons — after all, there are plenty of them 
around in a metal. We briefly consider key points on how to proceed in calculating 
the respective relaxation times for the three scattering scenarios. But first we men- 
tion a remarkable relation known for more than half a century prior to the devel- 
opment of quantum mechanics — Matthiessen’s Rule. This empirical rule, described 
by Matthiessen in 1862, 69 asserts that if several distinct scattering mechanisms are 
at play the overall resistivity of metals is simply the sum of the resistivities one 
would obtain if each scattering mechanism alone were present. For example, for 
two distinct scattering mechanisms — scattering on impurities and on lattice vibra- 
tions — the resistivity can be written as 



1 m* 

P = ~ = ^~ 
o ne~r 



ne~T(i) ne~T(2) 



P{\)+P(2)- 



(67) 



The essential point behind Matthiessen’s Rule is the independence of scattering 
mechanisms; i.e., the overall collision rate is the sum of the collision rates of the 
participating scattering processes. In the context of the relaxation time approxima- 
tion this immediately leads to a reciprocal addition of the relaxation times: 
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1 _ i i 

T T (1) T( 2 ) 



( 68 ) 



Extensive studies (for a review see Ref. 70) have been done to test the validity of 
Matthiessen’s Rule. Overall the rule seems to be reasonably robust, with deviations 
often less than 2%. However, there are situations where the rule breaks down, 
indicating that the scattering processes are not mutually independent after all. 
This happens when the relaxation time depends on the k vector; i.e., r = r(k). 
As discussed in Ref. 2, the resistivity is then proportional to a reciprocal average 
of the relaxation time 1 It, and Matthiessen's Rule implies 



1 _ i i 

r r (1) r (2 ) 



(69) 



On the other hand, averages performed on Eq. (67) require 

6 = (z L ) + (z L )- 

T T (l) T (2) 



(70) 



The formulas are not equivalent unless r(i) and t {2) are independent of k. Since 
there is a small but finite spread of electron wave vectors that can interact with 
phonons, and T( 2 ) are somewhat interdependent and this, in general, leads to an 
inequality rather than an equality between the total and two partial resistivities: 

P>P(l)+P(2)- (71) 



This inequality implies a positive deviation from Matthiessen’s Rule; i.e., if we 
write for the total resistivity p tot = Pi mp + p(T), the temperature-dependent resistiv- 
ity always increases with increasing impurity resistivity p„, qi . This old dogma has 
been challenged by measurements of Rowlands and Woods, 71 who observed nega- 
tive deviations from Matthiessen’s Rule in samples of Al, Ag, and Pd that were 
mechanically strained. Thus, it depends on how the resistivity p imp is being in- 
creased. If the sample is more “dirty” (by stuffing more impurity into it), the 
deviation will always be positive. If, on the other hand, p imp is increased by strain- 
ing the sample, then the deviation will be negative. The origin of the negative 
deviations is believed to be electron-dislocation scattering. In mechanically de- 
formed samples electron-dislocation scattering is dramatically enhanced and be- 
cause this interaction process favors small-angle scattering in contrast to large-angle 
scattering predominant in the case of electron-impurity scattering, this leads to a 
reduction of p(T) with increasing p im p . 

Although it is not perfectly obeyed, Matthiessen’s Rule gives us a hope of being 
able to discuss electron scattering in metals in terms of more-or-less independent 
processes and focus on important characteristics and signatures of each relevant 
scattering mechanism. Before we do that, let us consider the form of an electron 
distribution function as it is perturbed by a thermal gradient. 



3.3. Electronic Thermal Conductivity 

An expression for the electronic thermal conductivity is given in Eq. (52). If we do 
not neglect the term K\K~ X K\ with respect to the term K 2 , we can write the 
thermal conductivity as 
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= ^ [Ki - K,K-'K,} = i [tf 2 - eSK\T\ = ^ - eSK u (72) 



where we used Eq. (56), that defines the Seebeck coefficient S. Forming now a ratio 
k/ctT (i.e., defining the Lorenz function L). we obtain 

k % - eSK\ _ AT S'A'-'A'i tt 2 2 

txT = e 2 I\ 0 T ~ e 2 T 2 K 0 eT “T 



(tf) 



-S' 2 = T„ - S 2 



(73) 



Thus, the Lorenz function L differs from the Sommerfeld value L a by a term equal 
to the square of the Seebeck coefficient. In good metals at ambient temperatures the 
Seebeck coefficient is a few // V/K ; i.e., the term S 2 is ~ 10 11 V 2 /K 2 , which is three 
orders of magnitude smaller than L a and can be omitted in most situations. Never- 
theless, this correction term will be important when we discuss the high-temperature 
electron-electron contribution to the electronic thermal conductivity of metals. 

It is instructive to illustrate the electron distribution that arises as a result of the 
presence of a thermal gradient. In particular, by comparing it with the distribution 
created by an electric field, we will be able to highlight the difference between the 
two and understand why thermal and electrical processes differ at low temperatures. 

From Eq. (35), assuming only a thermal gradient is present, we can write for the 
distribution function of electrons 

/(k) = /°(k) + (-|g) [^)rv(k)(— VT) = f°(T + rv(k)(— VT)), (74) 



where the second equality follows from the use of Taylor’s theorem and the identity 

df° = E-^C ( df° 
dT T 



{ dE J 



(75) 



The distribution is shown in Fig. lb. We note the following important points: 
For E > ( and for v(k)- (— VT) >0 we obtain /( k) > /°(k) 
v(kH-VT) <0 /(k) < /°(k) 



(76) 



For E < ( and for v(k)- (— VT) >0 we obtain /( k) < /°(k) 
v(k)- (— VT) < 0 /( k) > /°(k) 



Therefore, for E > electrons going in the direction of -VT (i.e., down the ther- 
mal gradient) are hotter and tend to spread the distribution, while electrons going 
in the direction of VT (i.e., up the thermal gradient) are colder and sharpen the 
distribution. 

This is clearly a very different situation compared to that encountered in Fig. la, 
where the distribution was perturbed by an electric field with the result that the 
entire Fermi sphere shifted in response to the applied field. 

The nature of the thermal distribution in Fig. lb presents an alternative mechan- 
ism for the electron system to relax. While the electrons may achieve equilibrium by 
being scattered through large angles across the Fermi surface (the horizontal pro- 
cess effective in relaxing the electron population disturbed by an electric field in Fig. 
la), in Fig. lb electrons may relax simply by undergoing collisions in which they 




46 



Chap. 1.2 • THERMAL CONDUCTIVITY OF METALS 



change their energy by a small amount and move through the Fermi level, but with 
essentially no change in the wave vector. These types of processes — vertical pro- 
cesses — provide a very effective means of relaxing thermally driven distributions 
and govern thermal conductivity at low temperatures. At high-temperatures (usual- 
ly including a room temperature) electrons scatter through large angles, the average 
of cos 6 is zero, and a single collision may relax the distribution. Because scattering 
through large angles (horizontal processes) becomes less and less probable as the 
temperature decreases, electrons need many collisions before the direction of their 
velocity is changed substantially. However, if they can pick up or give up a bit of 
energy, no such limitation exists and they can relax readily. Thus, we are starting to 
appreciate that at low temperatures the relaxation time governing electrical con- 
ductivity will not be the same as the one governing thermal conductivity. The 
relaxation time in thermal processes will be considerably shorter because the elec- 
trons have an additional and very effective channel through which they can relax — 
the inelastic vertical process. This is the reason why the Wiedemann-Franz law 
breaks down at low temperatures — the electrical and thermal processes lack a 
common relaxation time. 



4. SCATTERING PROCESSES 

Electrons in a solid participate in three main interaction processes : they can scatter 
on lattice defects, they interact with lattice phonons, and they interact with each 
other. We have already noted that because of a small but finite spread of electron 
wave vectors, Matthiessen’s Rule is not perfectly obeyed and interference effects, 
especially between the impurity and electron-phonon terms, lead to deviations from 
the rule. Such a slight nonadditive nature of scattering processes arises from the 
energy dependence of the electron-phonon interaction. As one would expect, de- 
viations from Matthiessen’s Rule also affect the electronic thermal resistivity and, in 
fact, might be more pronounced there than in the case of electrical resistivity. The 
reason is that the electrical resistivity is more immune to small changes in the 
electron wave vector than is the electronic thermal conductivity. Nevertheless, the 
statement must be qualified in the following sense: (a) As pointed out by Berman, 1 
deviations from Matthiessen’s Rule are usually largest in the regime where the 
scattering rates for the electrical and thermal processes differ most — at intermediate 
to low temperatures. Since in this regime phonons are much less effective in limiting 
electrical resistivity than the electronic thermal resistivity, the fractional influence of 
the same amount of impurity will be larger in the former case, (b) Experimentally, it 
is far more difficult to design and implement high-precision measurements of ther- 
mal conductivity than to do so for electrical resistivity. Thus, while there is no 
shortage of the data in the literature concerning deviations from Matthiessen’s 
Rule in electrical resistivity, the corresponding data for the thermal conductivity 
are rather sparse. 

We now consider how the different interaction mechanisms give rise to electrical 
and thermal resistivities. 

4.1. Impurity Scattering 

Impurities and lattice imperfections are present to a greater or lesser extent in all 
real crystals, and they are an impediment to the flow of electrons. Scattering of 
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electrons on impurities is especially manifested at low temperatures, where the 
strength of the competing scattering mechanisms is weaker. In general, impurities 
are heavy objects on the scale of the electron mass, and, provided the impurity does 
not have an internal energy structure such that its energy levels would be closely 
spaced on the scale of kn r J\ the impurity will not be able to absorb or give energy to 
a colliding electron. Consequently, impurity scattering is considered a purely elastic 
event. Since there is no change in the energy of an electron between its initial and 
final states, the collision probability Wku is zero unless /: ( k ) = /:Tk'). If U is the 
interaction potential, first-order perturbation theory gives for the collision prob- 
ability : 

W w = ^6[E( k) - £(k')]|(k|[/|k')| 2 . (77) 

During linearization of the Boltzmann equation, we invoke the principle of detailed 
balance, which requires that Ww = Wwk- Thus, writing for the number of states of 
a particular spin in volume V, Vdk/( 2 tt) 3 , simplifies Eq. (32) to 

(78) 

The terms [1 -/(k)] and [1 -/(k')] fall out and the exclusion principle does not alfect 
the rate of change due to collisions. In order to relate the transition probability 
Wkk 1 to the relaxation time r, we must use a formula derived for the modified 
distribution function /(r,k,f). The simplest case is that corresponding to an electric 
field, and let us say we take it along the x-axis; i.e., 

df° 

f = f° + Teev x — . (79) 



From the definition of (df/dt) m n we have 




— e £ v x 



qr 

dE ‘ 



(80) 



Substituting for /(k) in Eq. (78), and equating the right-hand side with that of Eq. 
(80), we get for the relaxation time 

/ TT 

- = —J dk'Wiw [l -^1 f |(k|[/|k , )| 2 (l -cosd) sm0d9. (81) 

T (27 t) J { V x ) 27r?r J 



In deriving the second equality in Eq. (81), we use the scattering geometry in Fig. 2. 
Equation (81) indicates that the scattering frequency is weighted with a factor 11- 
cos 6). This means that not all scattering processes are equally effective in altering 
the direction of the electron velocity. Much more weight is given to those processes 
that scatter electrons through a large angle. 

Eet us consider the case of a static charged impurity of valence Z that occupies a 
lattice site. The conduction electrons “feel” the presence of this impurity through a 
screened Coulomb interaction. This arises because the conduction electrons will try 
to respond to the electric field of the impurity and will distribute themselves so as to 
cancel the electric field at large distances. In the simplest case, modeling the electron 
gas as a free-particle system, the screened potential has the form of the Yukawa 
potential : 
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U(r) 



Ze 2 

4ir£ 0 r 



-k T Fr 

C 7 



(82) 



where the parameter Utf is called the Thomas-Fermi screening parameter and its 
inverse (Jctf) is often referred to as the screening length. The Thomas-Fermi 
screening parameter is given solely by the electron density and can be written as 



k T F = 




(83) 



where fcp is the Fermi wave vector (or momentum) [fcp = (37 r 2 ™) 1 / 3 ] and a 0 = h 2 /me 2 
= 0.529xl0~ 10 m is the Bohr radius. For the electron densities in metals, the 
screening length is on the order of lattice spacing; i.e., electrons in metals are 
very effective in screening the electric held of impurity ions. The matrix element 
for this potential, <k'|[/ im? ,|k>, representing electron scattering from a state k to a 
state k' is the Fourier transform of the screened potential: 

, Ze 2 

k} = T. . , - (84) 



(k ' £/, 



imp | 



'°(l 



k-k'| 



rvrpF 



V 



Substituting this matrix element into Eq. (81), working out the integral, and in- 
troducing the density of impurities n, = N,/V, one obtains for the impurity resis- 
tivity , Pimp • 



m 2irm rii EpZ 2 ( 4 kip 

P' mp ne 2 T 3 ne 2 n h ’ 



(85) 




FIGURE 2 Scattering configuration in three dimensions. The electric field is assumed to be in the x- 
direction. Magnitudes of both initial and final electron wave vectors are Ap. A current change results 
from changes in the electron velocity with respect to the direction of the electric field. Assuming free 
electrons, the direction of the electron velocity coincides with the direction of its wave vector. Conse- 
quently, after averaging over <fi, the change in the current is given by (k x —k! x )lk x = 1 - v/Jv x = (1-cos 6). 
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where the function F is 

= 2[ln(l + x) -x/(l +x)] 
x 2 



( 86 ) 



It is of interest to estimate the impurity resistivity per 1% of impurity. As an 
example, we take copper where the Fermi energy is about 7 eV and the electron 
density is 8.5xl0 28 m -3 . Equations (83) (86) then yield the impurity resistivity 
contribution of 1.8x10 s m per 1% of monovalent (Z = 1) impurity. This result 
is in a reasonable agreement with experimental values (see Table 3), although only 
two entries (Ag and Au) can be considered monovalent solutes. 



The scattering cross section calculated with the Born approximation often tends 
to overestimate the cross section. The partial wave method is considered a better 
approach. In this scheme a charged impurity represents a spherically symmetric 
potential for which the Schrodinger equation has an asymptotic solution of the 
form 

ipk,i ~^sin ffcr-y + Tftj Pi(cos9). (87) 



Here P{( cos 6) is a spherical harmonic of order i. By matching the functions ipk,e to 
a plane wave (representing the initial state), one can write for the differential 
scattering cross section: 



a(9) 



1 

k 2 



Y (21 + \)e lr}l smr][Pi(cos9) . 
1=0 



( 88 ) 



TABLE 3 Resistivity increase in Cu per 1 at.% of Solute.” 



Impurity 


A pic 

(10-®S1 mil at.%) 


Impurity 


A pic 

(10 8 C2 m/1 at.%) 


Ag 


0.14 


Mn 


4.3 


A1 


1.25 


Ni 


1.14 


As 


6.5 


Pd 


0.84 


Au 


0.53 


Pt 


1.86 


B 


1.4 


Sb 


5.4 


Be 


0.64 


Se 


10.6 


Co 


5.8 


Si 


3.25 


Fe 


14.5 


Sn 


2.8 


Ge 


3.7 


V 


5.8 


In 


1.2 


Zn 


0.25 



“ Data are taken from extensive tables in Landolt-Borstein, New Series. 10 
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The phase shifts must satisfy the Friedel sum rule, 

Z = -J2(2l+l)m(k F ), (89) 

n i 

where Z is the valence difference between the impurity and the solvent metal. 

It is important to note that static lattice defects (such as charged impurities) 
scatter all electrons with equal effectiveness; consequently no temperature depen- 
dence arises in electrical resistivity due to static lattice defects. This contrasts with 
phonon transport, where the static defect scattering affects phonons of different 
frequencies differently, leading to a temperature dependence. 

We conclude this section by recapitulating the main points : Because the electrons 
participating in transport are highly degenerate and are contained within a narrow 
band of width ~/c s T around the Fermi level, as they encounter a static lattice defect 
they will scatter elastically. A static lattice defect “looks” the same to all electrons 
(a very different situation from that of phonons), and the relaxation time will be 
temperature independent. Of course, different static crystalline imperfections (either 
point imperfections, such as interstitial atoms and vacancies, or extended imperfec- 
tions, such as dislocations) scatter with different strength and effectiveness. Relaxa- 
tion times relevant to each specific situation have been worked out and can be 
found in the literature. 3,72 The essential point concerning all of them is that, in 
each case, the relaxation time is independent of temperature. Collisions with static 
imperfections result in a large-enough change in the electron wave vector so that 
both electrical and thermal processes are affected to the same degree. Electrical 
resistivity that is given purely by electron scattering on static imperfections is often 
called the residual impurity resistivity. It is the resistivity that a metal would attain 
at very low temperatures. The common relaxation time for electrical and thermal 
transport when electrons scatter on static imperfections means that the Wiede- 
mann-Franz law is valid. From this we glean that the impurity-dominated electro- 
nic thermal conductivity will be expected to follow linear temperature dependence 
and its inverse — thermal resistivity W imp — will be inversely proportional to tem- 
perature : 

HWT) = §r = %b (»> 



where Bimp — 

4.2. Electron-Phonon Scattering 

The most important process that limits the electrical and thermal currents has its 
origin in a particular kind of a crystal lattice imperfection — the disturbance of 
perfect periodicity created by a vibrating lattice ion. Just as an electron scatters 
when the lattice periodicity is interrupted by the presence of an impurity at a lattice 
site, the disturbance arising from vibrating ions (phonons in the language of quan- 
tum mechanics) leads to a deflection of an electron from its original path. The 
major difference between impurity scattering and electron-phonon scattering is 
that the former represents a purely elastic process, whereas the latter involves 
emission or absorption of a phonon, i.e., a small but finite change in the electron 
energy. From a practical perspective, one has some control over the impurity 
processes because one can, at least in principle, prepare samples of exceptional 
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crystalline quality with an exceedingly small amount of impurity in the structure. In 
contrast, there is nothing one can do about a perturbation created by a vibrating 
lattice — it is there to stay whether we like it or not. 

The topic of electron-phonon interaction is one of the pillars of solid-state 
physics. It presents a challenge that requires the full power of a quantum mechan- 
ical treatment of the problem. Rather than reproducing a lengthy and complex 
derivation that can be found in several texts, we will sketch the key steps and 
give results that are essential in the discussion of electrical and thermal transport 
in metals. Readers interested in a rigorous treatment of the electron-phonon inter- 
action may find it useful to consult, among many others, Refs. 3, 5, 63, and 73. 

In electron-phonon scattering we consider two distinct distributions, those of the 
electrons and phonons. We know how to treat each one separately-electrons are 
Bloch waves and obey Fermi-Dirac statistics, and vibrating metal ions show small 
displacements from their equilibrium positions and thus can be modeled as harmo- 
nic oscillators obeying Bose-Einstein statistics. In Chapter 1.1, in the discussion of 
heat transport by lattice vibrations, it was described how the coupled motions of 
atoms are transformed to form normal modes, each specified by its wave vector q 
and a frequency of vibration u q . These normal modes are quantized as harmonic 
oscillators, and the energy of a mode with frequency u q is (n + l/2)huj q , n being an 
integer. Excitations of the lattice are energy quanta hu> q , called phonons. The 
equilibrium distribution function, i.e., the mean number of thermally excited quan- 
ta hu q at a temperature T (the average number of phonons), is given by the Planck 
distribution function (Bose-Einstein statistics): 



N° = 

e hv q /k B T _ J 



(91) 



If the unit cell contains only one atom, there are only three acoustic modes 
(branches) of vibrations: one longitudinal with the atomic displacements parallel 
to the wave vector q, and two transverse, with the displacements perpendicular, to q 
and to each other. The dependence of the vibrational frequency uj q on the wave 
vector q, tv q (q), is called the dispersion relation. In the long-wavelength limit 
(q — >0), uj q =v s q, where v s is the speed of sound. If more than one atom resides 
in the unit cell, more modes arise and they are called optical branches (the number 
of optical branches is given by 3 (p - 1), where p is the number of atoms in the unit 
cell). Optical branches have non-zero frequency at q = 0. The group velocity of a 
mode, v g = cL> q /dq, indicates the speed with which wave packets travel in the crystal 
and transport thermal energy. Because the optical branches are rather flat, the 
derivative dto q /dq is quite small, which explains why optical branches are substan- 
tially ineffective as carriers of heat ; thus the task of transporting heat is assigned to 
acoustic phonon modes. Typical speeds of sound in metals are a few thousand 
meters per second compared to the speed of sound in air, 340 m/s. In describing 
the vibrational spectrum of solids, one often relies on the Debye model. A crystal is 
viewed as consisting of N lattice sites that give rise to 3 N acoustic vibrational 
modes. Moreover, the vibrational spectrum is truncated so that the highest fre- 
quency that can be excited — the Debye frequency u>d — is given by 



o 



D(u>)dio = 3 N. 



(92) 
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In Eq. (92) D(u>) is the phonon density of states equal to 

W = (93) 

where V 0 is the sample volume. The question is how to describe the interaction 
between the electrons and phonons. 

As pointed out, the feature that distinguishes the electron-phonon interaction 
from electron-impurity scattering is the quantum of energy hiv q that an electron 
absorbs or emits as it is scattered by a phonon. This is an important constraint 
because the electron’s energy is changed in the interaction, i.e., the scattering is 
inelastic. We assume that only one quantum of energy (one phonon) is involved in 
the scattering process, and therefore the scattering satisfies the requirement 

E(k) ±fiu = E(k!). (94) 

If the wave vectors of the incoming and outgoing Bloch states are k and k', and the 
phonon wave vector is q, the condition the wave vectors must satisfy (one can view 
it as a selection rule) is 

k ± q = k' + G (95) 

where G is a reciprocal lattice vector. Equation (95) delineates two distinct kinds of 
processes : interactions for which G = 0 are called normal or N-processes ; interac- 
tions where a nonzero vector G enters are known as umklapp or U-processes. In 
essence, Eq. (95) is a statement of the conservation law of momentum (p = hk); in 
the case of U-processes the vector k + q reaches over the edge of the first Brillouin 
zone, and the role of a reciprocal vector G is to bring it back to the first zone. As 
Fig. 3 illustrates, U-processes are exceptionally effective in generating resistance 
because they are capable of nearly reversing the direction of the electron wave 
vector. Of course, this necessitates a presence of phonons with a certain minimum 
wave vector which are more plentiful at high, rather than low, temperatures. There- 
fore, U-processes typically dominate at high-temperatures, and their influence 
weakens exponentially as the temperature decreases. 

Now, because the lattice ions are much heavier than electrons (M /m e ~ 10 3 ), it is 
assumed that the electrons quickly adjust to changes in the ionic positions. Con- 
sequently, within this concept, generally known as the Born-Oppenheimer approx- 
imation, the actual electron-phonon interaction is treated as a perturbation on the 
Bloch states that are calculated under the assumption that the ions are in equili- 
brium positions. 

Considering electron-phonon interaction, one might think that the lower-energy 
cutoff (some kind of minimum electrical resistivity) would be given by the zero- 
point motion of the lattice ions, i.e., by the term huj(q)/2 being emitted or absorbed 
by an electron. This, however, would be an incorrect viewpoint. To consider zero- 
point motion, we necessarily assume T = 0. But at absolute zero temperature 
neither phonon absorption nor phonon emission by electrons is possible. Absorp- 
tion does not happen because n q = 0, and phonon emission cannot take place 
because an electron would lose energy in such a process. Since at T = 0 the electron 
gas is perfectly degenerate, there is no available electron state below E F to accom- 
modate the electron. 
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FIGURE 3 (a) Electron-phonon N-process k+q=k'. (b) Electron-phonon U-process k+q=k'+G, where G 
is the reciprocal lattice vector. Note a minimum phonon wave vector q m that is necessary for the U- 
process to take place. At low temperatures such wave vectors may not be available. 



In order to describe the interaction between an electron and the lattice ions, let us 
assume that g(k,k',q) is a function that couples the Bloch states k and k' to the 
phonon of wave vector q via Eq. (95). We would like to find out the form of this 
coupling function, assuming a simple model that neglects any role of U-processes. 
This is the same simplification made by Bloch in the derivation of his Bloch- 
Griineisen formula for the electrical conductivity. Consider an electron situated 
at point r that “senses” a lattice ion located at point R, through a potential 
E/(r - Rj). The position of the ion R, can be specified with respect to its equilibrium 
position R“ and the displacement <5R, : . The potential energy of an electron in the 
field of all ions of the lattice is the interaction Hamiltonian H e - i on : 

U Q - io„ = 5] U [r - (R“ + SRi) ] , (96) 

i 

where the summation runs over all ion sites. The potential G[r - (R“+<5R ( ;)] has an 
electrostatic origin except that electrons have an excellent ability to screen the long- 
range tail of the Coulomb potential. One of the simplest forms of screening — 
Thomas-Fermi theory — assumes a slowly varying potential as a function of r given 
in Eq. (82) together with its screening parameter k T F hi Eq. (83). The required 
smoothness of the potential must be on the scale of a Fermi wavelength, which, 
stated equivalently, requires q <<fcf. i.e., the calculation using the Thomas-Fermi 
screening potential will be reliable only in the long-wavelength limit of lattice 
vibrations. Since we assume small displacements from the equilibrium position, 
we can expand the potential and, keeping only terms to first-order in <*>R,, we write 

K— ion = XI U(r - R° - 6Ri) = ]T U( r - R“) - SR, VC/(r - R") . (97) 

i i i 

The first term in the expansion is a constant and represents the potential to which 
the electrons are subjected when the ions are in their equilibrium lattice positions. 
This leads to the Bloch description of electrons in the crystalline lattice that we 
assume anyway. Thus, the electron-phonon interaction is described by the second 
term in Eq. (97). At this stage one can introduce the phonon coordinates Q qa via 
the displacements <SR* : 
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6R i; = TV'/ 2 ^ Q qa e qa e 

qa 



«q R" 



(98) 



where N is the number of lattice ions in a volume V 0 and e qa are unit polarization 
vectors chosen to be parallel or perpendicular to q when q is along a symmetry 
direction. 

From this point, one would usually resort to using the second quantization 
formalism and express the phonon coordinates in terms of creation and annihila- 
tion operators. This is an elegant approach and one that keeps the terms easily 
tractable. Because not all readers may be familiar with this technique, we shall not 
follow this road. We simply recognize that we must calculate matrix elements of the 
form <k'|<*>R,VE/(r-R°)|k> and sum them up for all ions i. The actual computation 
is given in Ref. 4. 

As noted, we only consider N-processes. If we further simplify the problem by 
treating electrons as free electrons (i.e., if we use a plane-wave exp(ikr)/(V 0 ) 1,/2 to 
describe an electron in state k), the coupling function for the electron-phonon 
interaction becomes a function of q = k/ - k only; thus 



• ? 
ie z 



--- / Nh 

q (1 + k^ F /q 2 )£ 0 q 2 V 0 V 2Mut 






(99) 



Here, z is the charge state of an ion, and the subscript a is dropped because the 
polarization is assumed to be either longitudinal or transverse to the wave vector q. 
We need the square of the modulus of C, which in the long-wavelength limit is 



, 1 2 2 n 2 zEphcjq 

'9 1 9-0= 3 v 2 k ^ Vo M ' 



( 100 ) 



Eq. (100) is usually written in the form 



I f< I 2 _ \ ^9 

| 9 g— >0 A 2N(0)V o ' 



( 101 ) 



where 



N( 0) 



mkF 
2n 2 Tr ’ 



( 102 ) 



is the density of states per spin at the Fermi level, and the dimensionless parameter 
A is 

\- 2z E F m e E F ( gjQ 2 ( c ) 2ffl e 

3 (fc B 0 ,) 2 M ~ 3 ( k B e D f UfJ UJ m ' 

The temperature 0 S in Eq. (103) is related to the longitudinal sound velocity v s via 
the equation 

0 S = v s hk F /k B - (104) 

The temperature 0 S defined in Eq. (104) and the Debye temperature 0d, in terms of 
which most of the asymptotic forms of the transport formuli and their integration 
limits are usually stated in the texts, are related by 
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TABLE 4 Temperatures 0 S and Bo Together with the Parameter A of Eq. (103) for Selected Metals. a 



Metal 


0« (K) 


6d (K) 


Bp (K) 


A (dimensionless) 


Na 


220 


150 


200 


0.47 


K 


150 


100 


114 


0.25 


Cu 


490 


315 


330 


0.16 


Ag 


340 


215 


220 


0.12 


A1 


910 


394 


395 


0.90 



Q Values of O s and A are from Ref. 4, those of &o are from Ref. 74, and 6r is from Ref. 75. 



0s _ Vs_kF 
do c qo 



(105) 



Here c is the speed of sound in the Debye model, which is understood to be the 
average of the long-wavelength phase velocities of all three acoustic modes, not just 
the longitudinal one v s . Since for free electrons Uf - (2/Z) 1/3 q D , where Z is the 
nominal valence and q/, is the Debye wave vector, 0 S could be larger (perhaps as 
much as a factor of 2, depending on Z) than the Debye temperature d D (see Table 4.). 



Summarizing the approximations made to obtain the result in Eq. (101), we 
have: 



a. Interaction processes were restricted to involve only N -processes 

b. Thomas-Fermi screening theory was used 

c. Electrons were treated as free electrons. 



Such simplifications, of course, have consequences in terms of the validity of the 
calculation. Even in the simplest of metals (alkali and noble metals), U-processes 
cannot be completely neglected. Moreover, the coupling constant is appropriate for 
the long-wavelength limit only and not for phonons of higher frequency. Finally, 
Fermi surfaces of metals are, in general, quite complicated and electrons are not 
exactly free electrons. In spite of these limitations, the form of the electron-phonon 
coupling constant derived is used rather indiscriminately in most calculations. The 
reason is obvious — such a treatment yields a rather simple and compact result, and 
alternative approaches lead to very messy calculations, if possible at all. 

With the coupling constant given by Eq. (101), one may proceed to calculate the 
collision integral (df(k)ldt) co \\ and inquire about the form of the relaxation time 
describing the electron-phonon interaction: 

[ %) = "irE l c 'l 2 { [Mi ~ Mi + N -i) - M 1 - - E k + fo- q ) 

\ Ol ) coll n 

+ [fk( 1 - h)N q - Ml - /*) (1 + N q )]6(E„ -E k - huj q )}. (106) 

The factors such as //,■( I - f)/) enter because of the Pauli exclusion principle, and N q 
and N_ q stand for phonon distributions under the processes of phonon emission 
and absorption, respectively. Of course, when both electrons and phonons are in 
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equilibrium, the collision term must vanish. When the distributions are disturbed by 
an electric field or thermal gradient, the deviations from the equilibrium population 
can be expressed with the aid of the deviation functions ipk for electrons and <1*,, for 
phonons : 



fk = /° + /°(l -/“)#, 


(107) 


N q = N° + N°(l + N°)<j) q . 


(108) 



Computations are considerably simplified if one makes the following approxima- 
tion. Provided the phonon-phonon and phonon-impurity scattering processes are 
very frequent so that they relax the population of phonons effectively, the electron 
will interact with a phonon that is very close to equilibrium. Hence, one sets the 
deviation function of phonons to zero. With the identity equation that relates 
equilibrium distribution functions of electrons and phonons, 

f°m 1 - f(E+hw)]N° M = f°{E + hw)[\ - f°(E)](l + N°(u>)), (109) 



one can linearize Eq. (106). Introducing a relaxation time 7>., one simplifies Eq. 
(106) to 
1 2t r. 



Tk 



= jrJ2 l^l 2 S ( E k - E « ~ H) [ N ° q + 1 - ft) + K E k - E « + foq) ( N ° q + ft 



( 110 ) 



The coupling function depends only on the difference between the wave vectors of 
the outgoing and incoming electrons, q = k' - k. Assuming again that the electrons 
are described by a free-electron model and converting from summation to integra- 
tion, we obtain for the high- and low-temperature limits the following temperature 
dependences for the electron-phonon relaxation rate: 

1 It, cx T for T > 9 , 

' k , ~ ° 111 
cxT for T « 6 d . 



The high-temperature result — the relaxation rate being a linear function of tem- 
perature — could have been anticipated based on a classical argument utilizing equi- 
partition of energy, which is valid since we consider high-temperatures. One would 
expect the resistance to be proportional to the mean square displacement of an ion 
in a given direction. An ion of mass M with acceleration d 2 xldt 2 is subjected to a 
restoring force -bx so as to satisfy the equation 

Mx + bx = 0. (112) 

Since b/M = oj 2 = 4ir 2 f 2 , where / is the frequency, we can write for the mean 
potential energy of this oscillator at high-temperatures : 



The Debye temperature qD is the usual upper limit of the summation in Eq. (110) and. thus, the 
temperature delineating the asymptotic regions of transport. However, this is applicable only for highly 
degenerate systems such as metals. For semimetals a more appropriate criterion is q / 2 < Hf- Moreover, if 
the model assumes that electrons couple only to longitudinal phonons, a somewhat higher temperature, 
such as ©,, in Eq. (104) might be more appropriate. 
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i txl='-k B T. 

The resistance then becomes 

- k B T k B T h 2 T T 

(XX b 4ir 2 f 2 M 4n 2 Mk B 0 D ^ M0 2 D 



(113) 

(114) 



In order to appreciate what happens to an electron when it interacts with a 
phonon at high-temperatures, it is instructive to consider the largest possible 
change in the direction of an electron wave vector in such an electron-phonon 
encounter. If we assume a free-electron-like monovalent metal (one free electron 
per lattice ion) consisting of N atoms, the Brillouin zone must accommodate N 
distinct ^-values (modes of vibrations), while the Fermi sphere must be large en- 
ough to fit (M2) ^-states (M2 rather than N because each state can have two 
electrons, one with spin up and one with spin down). Thus, the volume of the 
Brillouin zone will be twice as large as the volume of the Fermi sphere, which 
implies that the ratio qD/k F = 2 1//3 = 1.26. Therefore, the maximum change in 
the electron wave vector is Ak = q B = l.26k F . From Fig. 4 it follows that the 
largest possible angle 9 is given by sin(6V2) = (q/2 )/k F . Inverting this expression, 
we obtain the maximum scattering angle d max = 2 sin _1 (q/2/cF) = 78°. Of course, 
this is an idealized case not fully applicable even in alkali metals. Nevertheless, the 
estimate shows clearly that high-energy phonons (which are plentiful at high-tem- 
peratures) can, in a single collision, produce a large change in the direction of an 
electron with very little change in its energy. Because high-temperature phonons can 
relax out-of-equilibrium populations of electrons created by either an electric field 
or a thermal gradient in a single scattering event, there is one common relaxation 
time for both electrical and thermal processes. 

As for the low-temperature limit of Eq. (Ill), it states that for temperatures well 




FIGURE 4 Relationship between the Fermi and Debye spheres for free electrons: q B = 2'^k F = 1.26kf. 
The maximum scattering angle d masL is given by 9 Iuax = 2 sin^ 1 (q B l2kp) = 78°. 
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below 0 D the electron-phonon scattering rate declines as T 3 . This must not be 
confused with the rate at which the current density diminishes because the two 
rates are not the same. What comes into play is the efficiency of scattering. It 
expresses the fact that, as the temperature decreases, phonon wave vectors become 
smaller and therefore the scattering angle decreases. Thus, more collisions are 
needed to make a substantial change in the direction of the electron velocity. We 
have seen in Eq. (63) that this is taken care of by introducing the efficiency factor 
(1 - cos 9). From Fig. 4 it follows that 



1 — cos 9 = 2 sin" - = 






(115) 



Because q ~ k B T/hv s , the factor (1 - cos 9) introduces a temperature factor T 2 . 
Thus, at low temperatures, the electrical resistivity is proportional to T 5 . This result 
can be derived rigorously by using Kohler’s variational technique to actually solve 
the Boltzmann equation. We will not reproduce the mathematical derivation, but 
merely state the final results. 



1 



9ttH 2 C 2 



Bv / 2( m* ) l/2 ea 2 Mk B 9 D E F y 2 



M)- Ms) ’*(?)■ 



(116) 



where 



and 



9Trh 2 {m*) l/2 C 2 
8 V2ne 2 a i Mk B 9 D E F y 2 



(?) 



9/T 



Jn = = 



( e x - 1 )“ 



~dx. 



(117) 



(118) 



Inverting Eq. (61) and substituting for the relaxation time from Eq. (116), one 
obtains the famous Bloch-Griineisen formula for the electrical resistivity: 76,77 



Pe-p{T ) 



m 



= A 






5 m ?)- 



(119) 



We indicate by subscript e-p that the resistivity arises as a result of the electron- 
phonon interaction. In most texts such a resistivity is referred to as the ideal 
resistivity — ideal in the sense that it would be the resistivity in the absence of 
impurities (i.e., the resistivity of an ideal solid). Note that different authors define 
the coefficient A differently. Apart from the form of Eq. (117), one often conies 
across a coefficient that is four times smaller, namely, A' = A/4, in that case, Eq. 
(119) is written with a prefactor of 4 A'. 

It is of interest to look at the asymptotic solutions for low and high-tempera- 
tures : 



Pe-p{T) = A 




\9 d /T) a 

4 



A 

4 



(a 



for T > 9 b 



( 120 ) 



and 
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p e _ p (T) = 124.4 A for T « 0 D . (121) 

In Eq. (120) we used the approximation J n (0/T) = ((VT)^ -1 ^™ -1 ^, valid for high- 
temperatures, and in Eq. (121) we substituted J 5 (o o) = 124.4, applicable at 
T « 6d ■ Thus, at high-temperatures the resistivity is a linear function of tempera- 
ture, while at low temperatures it should fall very rapidly following a T 5 power law. 
Note that the T 5 behavior of the resistivity is rarely observed. Apart from the 
obvious experimental challenges to measure resistivity at low temperatures with a 
precision high enough to ascertain temperature variations against a now large 
impurity contribution and a contribution due to electron-electron scattering (see 
the next section), we should keep in mind the simplifications introduced in the 
process of calculations, the main ones being spherical Fermi surfaces, N-processes 
only, Thomas-Fermi approximation, equilibrium phonon distribution. All of these 
do affect the final result. It is remarkable that with basically the same approach the 
theorists are able to calculate an electrical resistivity of even transition metals that 
agrees with the experimental data to within a factor of 2 (see Table 5). 



We note that by fitting the experimental resistivity data to the Bloch-Griineisen 
formula [Eq. (119)], we can obtain a characteristic temperature 0 R that plays a 
similar role as d D or the characteristic temperature associated with long- wavelength 
phonons, 0 S . The temperature d R is usually much closer to d D (see Table 4) than to 
0, s , indicating that transverse phonons do play a role in electron scattering. 

We now turn our attention to a relaxation time relevant to electronic thermal 
conductivity. This is obtained by performing similar variational calculations to 
those that led to Eq. (116) except that the trial function is now taken to be propor- 
tional to E — E f : 

J__ 9tt h 2 C 2 ( T\ 5 j ( x 

t k - 2\/l (m* ) 1/2 a 3 M k B 0 D E F ^l 2 U/J 5 ItJ * 



' 3 rfcf'i 2 (6_d\ 2 i Me/ry 

, ^ 2 Ud) ItJ 2 n 2 J 5 (0/T), 



( 122 ) 



TABLE 5 Comparison of Calculated and Experimental Room Temperature Electrical Resistivities of 
Several Pure Metals." 



Metal 


Ag 


Au 


Cu 


Fe 


Mo 


Nb 


Ni 


Pd 


Ta 


W 


Experiment 


1.61 


2.2 


1.7 


9.8 


5.3 


14.5 


7.0 


10.5 


13.1 


5.3 


Calculation 


1.5 


1.9 


1.9 


11.4 


6.1 


18.4 


5.2 


5.2 


13.2 


7.2 



“ Calculated values are from Refs. 78 and 79, and the experimental data are from Ref. 80. Resistivities 
are in units of 10~ 8 f2 m. 
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Electronic thermal resistivity arising from electron-phonon interaction (called ideal 
thermal resistivity) is then 



A 

~Ljf 




1 



1 Ji(0/T) 

2t r 2 J 5 (0/T) 



(123) 



where C v = (n 2 nk 2 B /2EF)T is the electronic specific heat. The ideal thermal resis- 
tivity consists of three distinct contributions. The first contribution is recognized as 
having the same functional form as the one in Eq. (119) for electrical conductivity. 
This term arises from large-angle scattering and, therefore, satisfies the Wiede- 
mann-Franz law. The second term is due to inelastic small-angle scattering (vertical 
processes) and has no counterpart in the electrical resistivity. The third term is a 
correction that accounts for situations where large-angle scattering can reverse the 
electron direction without actually assisting to restore the distribution back to 
equilibrium. This is due to the nature of the thermally perturbed distribution func- 
tion, where an electron can be scattered through large angles between regions 
having similar deviations (for instance, regions marked + in Fig. lb). 

Frequently, asymptotic forms of Eq. (123) are required for low and high-tem- 
peratures, and they are 






(124) 






' 0 D ' 



p e — P (L.T.) 3 (k F s 



L n T 7T“ 



m f^l =37.8^- m for T« 

Ud) ItJ L 0 9 D {q D ) { e D ) 



0d. (125) 



Equation (124) is just the first term of Eq. (123) in the limit of high-temperatures. 
The second and third terms in the square brackets of Eq. (123) are both very small 
because at high-temperatures J7/J5 — +(l/6)(0/T) 6 /(l/4)(#/7) 4 — >(2/3)(d/T) 2 — >0. The 
second equality in this equation follows from using Eq. (120). Thus, at high-tem- 
peratures, we expect the ideal electronic thermal resistivity to be temperature in- 
dependent and, because the relaxation times for the electrical and thermal processes 
are essentially identical, the Wiedemann-Franz law is obeyed, as Eq. (124) demon- 
strates. 

Equation (125) was obtained by considering only the second term in the bracket 
of Eq. (123) because this term rises rapidly as temperature decreases, while the first 
and third terms are constant (1 and 5082/124.4, respectively). The second equality 
in Eq. (125) is obtained with the aid of Eq. (121). The third equality follows from 
approximating J 5 (9d/T) = J 5 (oo) = 124.4. We note an important feature — the 
thermal resistivity at low temperatures follows a quadratic function of temperature 
and its inverse, the thermal conductivity n e - p (T), should be proportional to T~ 2 . 
This is indeed what is most frequently observed. Considering all the approximations 
made, the numerical factors should not be expected to be very precise, perhaps to 
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within a factor of 2 or 3. We also observe from Eq. (125) that, at low temperatures, 
the Wiedemann-Franz law is clearly not obeyed. 

Writing for the Lorenz ratio L=n!aT and substituting from Eqs. (119) and 
(123), we obtain 
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\6 d J 
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Asymptotic forms of the ratio L/L a at high and low temperatures are 



— = 1 for T > 0 D 
Ln 



(126) 



(127) 



<■»> 

Equations (127) and (128) state an important result. At high-temperatures, the 
Lorenz ratio attains its Sommerfeld value L 0 . In other words, the relaxation times 
for electrical conductivity and thermal conductivity are identical, and the Wiede- 
mann-Franz law holds. At these temperatures there are plenty of large-momentum 
(wave vector) phonons, and they scatter electrons elastically through large angles 
regardless of whether electron populations are perturbed by an electric field or a 
thermal gradient. At low temperatures, a very different scenario emerges. Elere the 
Lorenz ratio rapidly decreases with decreasing temperatures and the Wiedemann- 
Franz law is obviously not obeyed. Phonon wave vectors decrease with temperature 
and they are incapable of scattering electrons through large angles. Thus, many 
collisions are necessary to change the velocity of electrons and therefore relax the 
electron distribution brought out of equilibrium by an electric field. In contrast, it is 
much easier to relax a thermally driven disturbance in the electron population. A 
very efficient route opens whereby electrons just slightly change their energy (in- 
elastic processes) and pass through the thermal layer of width about the 

Fermi level, and the distribution is relaxed. This means that the relaxation time 
for thermal processes is considerably shorter than the one for electrical processes. 
The rate with which the two relaxation times depart as a function of temperature is 
given by Eq. (126), and a plot of L/L 0 as a function of reduced temperature is 
shown in Fig. 1.10 in Sec. 6. At low temperatures the ratio attains a quadratic 
function of temperature. The decrease in the Lorenz ratio is particularly large in 
very pure metals because the temperature range where Eq. (128) is applicable ex- 
tends to quite low temperatures before any influence of impurities is detected. 
Eventually, at some very low temperature, impurity scattering takes over and, since 
such processes are essentially elastic, the Lorenz ratio is restored. 



4.3. Electron-Electron Scattering 

Because of their high electron density, one might naively think that the electron 
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electron interaction is a very important component of the electrical and thermal 
conductivities of metals. Just a rough estimate of the Coulomb energy [e 2 /47re 0 r] of 
two electrons separated by an interatomic distance comes out to be about 10 eV — 
an enormous amount of energy that exceeds even the Fermi energy. However, this 
is a very misleading and incorrect estimate. The point is that electrons have an 
exceptionally good ability to screen charged impurities. We have noted that, in 
metals, the screening length is typically on the order of an interatomic spacing. 
This fact, together with the constraints imposed by the Pauli exclusion principle, 
results in quite a surprising ineffectiveness of the electron-electron processes in 
metals, especially as far as electrical transport is concerned. 

Interactions among electrons involve four electron states: two initial states that 
undergo scattering and two states into which the electrons are scattered. It should 
be clear that if only N -processes (processes in which the total electron momentum is 
conserved) are involved, the electron-electron interaction on its own would not lead 
to electrical resistivity because there is no change in the total momentum.* Thus, a 
possibility for a resistive process in interactions among electrons is tied to the 
presence of U-processes. But such processes are not too frequent because of the 
rather stringent conditions imposed by the Pauli exclusion principle. This is readily 
understood from an argument given in Ref. 2. Here we consider a metal at 7=0 
with fully occupied states up to the Fermi level and a single electron excited just 
above the Fermi level. Let us label this electron 1 ; its energy is E\ > E F . In order 
for this electron to interact, it must find a partner from among the electron states 
lying just below E F because only these states are occupied. We label one such 
electron 2, and its energy is E 2 < E F . The Pauli principle demands that after the 
interaction the two scattered electrons must go to unoccupied states and these are 
only available above the Fermi level; i.e., E[ > E F and E' 2 > E F . Energy in this 
scattering process must be conserved, 

E l +E 2 = E\+E' 2 , (129) 

and the wave vectors must satisfy 

k, + k 2 = k' 1 + k; + G. (130) 

If Ei is exactly equal to E F , then E 2 , E\, and E 2 must also equal E F . Thus, all four 
states would have to fall on the Fermi surface occupying zero volume in k-space 
and therefore yielding infinite relaxation time at 7=0 K. With E\ slightly larger 
than E F (see Fig. 5), there is a very thin shell of k-space of thickness \Ei—E f \ 
around the Fermi level available for final states E\ and E 2 . Because only two states 
(e.g., E 2 and E\) are independent rather than three [E 2 is fixed by Eq. (129) once 
Ei, E 2 , and E\ are chosen], the scattering rate is proportional to (Ei - E F ) 2 . 

If we now assume a finite temperature, the electron distribution will be slightly 
smeared out (by an amount of thermal energy k B T), and partially occupied levels 
will emerge in a shell of width k F T around E F . Since each one of the two inde- 
pendent states now has this enlarged range of width k F T of possible values avail- 
able, the scattering cross section is reduced by ( k F T /E F ) 2 ~ 1(R 4 at ambient tem- 
perature. Thus, at room temperature, the contribution of the electron-electron 
processes to the transport effects in bulk metals is negligible in comparison to other 



This is strictly correct only for isotropic systems. Any source of anisotropy will give rise to a resistivity 
contribution associated with N-processes. 
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FIGURE 5 Schematic of a normal electron-electron interaction. Electrons 1 and 2 with initial wave 
vectors ki and ki scatter from one another to final states with wave vectors ki' and ki'. The scattered 
electrons must go into unoccupied states (i.e., states above the Fermi energy). 



scattering mechanisms. If one takes into account static screening, the collision 
frequency for a purely isotropic electron-electron scattering rate has been 
evaluated 81 and is 



1 

Te-e 




(131) 



where 



l(x) =^(V£ 0 k 2 F /e 2 Y 



cos(d/2) 



x 

4 



tan 1 x + 



V2- 



itan 1 x\J 2 + x 2 



(132) 



The angular brackets denote an average over a solid angle, which is evaluated 
explicitly in Ref. 82. 

The important point for our purposes is the T 2 dependence of the scattering rate. 
As the temperature decreases, the collision rate diminishes. The T 1 variation applies 
to bulk, three-dimensional metals. If a metallic structure has lower effective dimen- 
sionality, a different temperature dependence results. For instance, for a one-dimen- 
sional wire the electron-electron scattering rate is a linear function of temperature 
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TABLE 6 Calculated Values of the Coefficient A ee (in units of 10 15 0 m K 2 ) and the Effectivenes 
Parameter A (percent), Together with the Available Experimental values of A ee . a 



Metal 


r Ja„ 


A ee (calc) 


A 


A ee (calc)/ A A ee (exp) 


Reference b 


Li 


3.25 


2.1 


0.054 


40 


30 


94 


Na 


3.93 


1.4 


0.035 


40 


1.9 


95 


K 


4.86 


1.7 


0.021 


80 


0.55^1.0 


95-98 


Rb 


5.2 


3.5 


0.028 


130 






Cs 


5.62 


6.7 


0.028 


240 






Cu 


2.67 


1.2 


0.35 


3 


0.28-1.21 


99-102 


Ag 


3.02 


1.1 


0.35 


3 


0.35T0.15 


100, 103 


Au 


3.01 


1.9 


0.35 


5 


0.5T0.2 


101 


A1 


2.07 


4.1 


0.4 


10 


2.8T0.2 


104 


Pb 


2.3 


16 


0.4 


40 






Fe 


2.12 








310 


105 


Co 










98 


106 


Ni 










340 


107 


Nb 


3.07 








23 


108 


Mo 










12.6 


109 


Ru 










27 


110 


Pd 










350 


111 


W 










4. 8-6. 4 


112 


Re 










45 


113 


Os 










20 


113 


Pt 










140 


114 


Nb 3 Sn 










70,000 


115 


Bi 


2.25 


80.000 






135,000 


116 


Graphite 










50,000 


117 



a Also included is the electronic density normalized to the Bohr radius, r s /a 0 . 

h Values of A for alkali metals are from Ref. 91; those for noble metals are from Ref. 92; those for the 
polyvalent metals are from Ref. 82. The calculated values A ee (calc) include a phonon-mediated correc- 
tion evaluated in Ref. 91. The theoretical estimate of A ee for Bi and graphite is based on Ref. 93. 
References are provided for the experimental values of A ee . 



while for a two-dimensional thin metal film the electron-electron interaction yields 
a contribution proportional to T 2 \n(T F /T). Since our ultimate interest is the elec- 
tronic thermal conductivity, and heat conduction measurements on lower-dimen- 
sional systems are exceedingly difficult to carry out, we limit our discussion to bulk 
metallic systems only. 

Equation (132) implies a reduction with decreasing temperature of the already 
weak electron-electron scattering rate. This, however, does not mean that we can 
completely neglect the influence of electron-electron processes. At very low tem- 
peratures, where the electron-phonon processes rapidly disappear (p e - P oc T 5 ), and 
in very pure metals, where the impurity resistivity contribution is small, the rela- 
tively slow T 2 dependence may propel the electron-electron processes into a posi- 
tion of the dominant scattering mechanism. 

4.3.1. Effect of e-e Processes on Electrical Resistivity 

The earliest experimental evidence for a T 2 term in the resistivity due to electron- 
electron interaction is in the data on platinum measured in 1934 by de Haas and de 
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Boer. 83 Shortly afterward, Landau and Pomeranchuk 84 and Baber 85 drew attention 
to the role of electron-electron scattering in the transport properties of metals. 
Prior to World War II, it was well established that, among the metals, transition 
elements are the ones that display an identifiable and measurable T 2 term in their 
low-temperature resistivity. A textbook explanation why transition metals behave 
this way was provided by Mott 86 in his model based on the sd scattering process. A 
characteristic feature of transition metals is an incomplete d band with high density 
of states, in addition to partially filled s and p bands of more or less ordinary (i.e., 
mobile) electrons. As the s electrons try to respond to an external field, they scatter 
preferentially into the high-density d states and acquire the character of the “slug- 
gish” d electrons. This leads to a rather dramatic decrease in the electric current 
density. Electron-electron scattering of this kind does not even require U-pro- 
cesses; the ordinary N -processes perfectly suffice to generate the large resistivity 
observed in the transition metals (see Table 1). Because more recent studies indicate 
that the density of states in transition metals is not that much larger in comparison 
to the density of states in noble metals, an alternative explanation has been put 
forward 87 to explain the large T 2 term in the resistivity of transition metals. It 
stresses the importance of a realistic shape for the energy surfaces and the use of 
the t matrix rather than the Born approximation in calculations of the electron- 
electron scattering cross section. Nevertheless, the intuitive appeal of Mott’s picture 
is hard to deny. 

In the mid-1970s, following the development of high-sensitivity, SQUID-based 
detection systems that allowed for an unprecedented voltage resolution in resistivity 
measurements, major experimental efforts confirmed that the electron-electron in- 
teraction has its unmistakable T 2 imprint on the low-temperature resistivity of most 
metals (Table 6). Moreover, the interest in electron-electron processes has been 
hastened by the development of localization theories, and throughout the 1980s 
much has been learned about electron-electron interaction in lower-dimensional 
structures. Readers interested in these topics are referred to Kaveh and Wiser. 88 

Among the outcomes of the intense theoretical effort were two new perspectives 
on electron-electron scattering in metals. The first concerns situations where a 
metal is exceptionally impure, i.e., when the elastic (impurity) mean free path is 
reduced to a size comparable to the interatomic spacing. In this regime of transport, 
the usual Fermi liquid theory that is at the core of all our arguments breaks down 
and must be replaced by a theory that takes into account the effect of disorder on 
the electron-electron interaction. 89,90 Specifically, under the conditions of strong 
disorder, in addition to the Fermi liquid term T 1 of Eq. (131), there appears a 
new term 

— = C(£ e )(k B T) d ' 2 , (133) 

'T'e—e 

where C depends on the mean free path l e and d is the dimensionality of the 
metallic system. Thus, the overall scattering rate due to electron-electron processes 
in very “dirty” metals is 

= AeT 2 + C(Q(k B Tf 2 = A ee T 2 + i ( J-) V2 3/ V 2 , (134) 

where in the second equality the constant C(l e ) is written explicitly. For reasonably 
pure metals ( k F l e >>1) the second term in Eq. (134) is negligible in comparison to 
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the first term. However, for a very impure metal for which the product kyl,. ap- 
proaches the Yoffe-Regel criterion, kyl, ~1 (near a metal-insulator transition), the 
second term in Eq. (134) dominates. We shall not consider such highly disordered 
metals. Rather, we focus our attention on systems that can be treated in the spirit 
of the Fermi liquid theory — essentially all bulk metals. Therefore, we take it for 
granted that the scattering rate for electron-electron processes is simply just the 
first term in Eq. (134). 

The second important feature emerged from a rigorous many-body calculation of 
MacDonald 118 and relates to a strong enhancement in the electron-electron colli- 
sion rate observed in the low-temperature electrical resistivity of metals, which 
ultimately undergo a transition to a superconducting state. This enhancement arises 
from the phonon-exchange term in the “effective’' electron-electron interaction. 
The reader may recall that in a conventional picture of superconductivity (the 
BCS theory) the mechanism of electron pairing is an attractive interaction between 
electrons via an exchange of virtual phonons. In nonsuperconducting metals (e.g., 
alkali metals and noble metals) this attractive interaction between electrons is 
smaller than their Coulomb repulsion. In metals that do superconduct, the attrac- 
tive interaction is the dominant interaction. This effect leads to a spectacular en- 
hancement (one to two orders of magnitude) in the electron-electron scattering rate 
of the polyvalent metals such as aluminum or lead. However, at high-temperatures 
this enhancement is substantially washed out. The limited range of temperatures 
where the effect shows up implies that one should be cautious when comparing the 
influence of electron-electron processes on the electrical resistivity to its influence 
on the electronic thermal conductivity with the latter accessible only at high-tem- 
peratures, as we discuss shortly. 

With the electron-electron scattering term given by Eq. (131), the electrical re- 
sistivity Pe—e and the electronic thermal resistivity W e - e (the subscript e-e refers to 
the electron-electron processes) can be written as 



Pe-e(T ) = A ee T 2 


(135) 


and 




We-e = B ee T. 


(136) 


The form of Eq. (136) is consistent with the Wiedemann-Franz law, in which case 
the ratio A ee /S ee plays the role of L. There are indications 119 that the value of L for 
electron-electron processes is close to 1 .1 x 1 0 8 V 2 /K 2 . 

An experimenter faces a daunting task when measuring the coefficients A ee and 
B ee ■ Writing for the overall electrical and thermal resistivities 


P — Pimp + Pe—e(T) + p e — p (T ), 


(137) 


W = Wimp(T) + W e - e (T) + W e - p (T), 


(138) 



one must be able to resolve the electron-electron terms against the contributions 
arising from impurity scattering and from the electron-phonon interaction. Let us 



Depending on what impurities are present in a given metal, there might be additional terms in Eq. 
(137) representing Kondo-like resistivity and possibly a T-dependent electron-impurity term, which, 
unlike in Eq. (85), arises due to inelastic scattering. 
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illustrate how challenging the task is for the easier of the two cases — determination 
of the coefficient A ee . 

In copper, for example, the impurity resistivity in the purest samples of copper is 
on the order of 10 12 <1 m. Electrical resistivity due to electron-electron processes is 
on the order of 10 16 S7 m at 1 K. Even if one measures the electrical resistivity at 
very low temperatures (T < 1 K) in order to eliminate a contribution from the 
electron-phonon scattering, the ratio p e - e ( 1 R)/p™p ~ 10~ 4 . To resolve this with an 
accuracy of 1%, one is called upon to make precision measurements at the level of 
one part per million. Such precision was impossible to achieve prior to the devel- 
opment of high-sensitivity SQUID-based voltage detectors (e.g., Ref. 120). How- 
ever, once these supersensitive devices became available and were incorporated in 
dilution refrigerators or helium-3 cryostats, a wealth of data appeared in the lit- 
erature ready to be analyzed. 

Even in the purest of specimens, one eventually enters a temperature regime 
where p imp becomes the dominant resistive process. Because of its presumed tem- 
perature independence, one can eliminate the impurity contribution by simply tak- 
ing a derivative of Eq. (137) with respect to temperature and isolate the electron- 
electron term by writing 

1 dp = 1 dp imp 1 dp e - e (T) 1 dp e _ p (T) = 1 dp e _ p (T) 

2 TdT 2 T dT 2 T dT 2 T dT ee 2 T dT ' 1 ’ 

If one now extends measurements to sufficiently low temperatures to suppress the 
electron-phonon term, one obtains A ee . This scheme works well regardless of 
whether the electron-phonon term is a power law of temperature or, as in the 
case of alkali metals, an exponentially decreasing contribution on account of a 
phonon drag effect. 121 

The values of A ee for transition metals (Table 6) are by far the largest, and thus 
the data have been collected with relative ease for most transition elements. Inter- 
pretation of the data, however, is often complicated by the presence of other 
interaction mechanisms, most notably electron-magnon scattering which also man- 
ifests its presence by a T 2 temperature dependence. It is interesting to compare the 
coefficients A ee (experimental and theoretical values) for alkali metals with those for 
noble metals. The data in Table 6 indicate that there is little difference between the 
two, at least if one focuses on the lighter alkalis. This is a rather fortuitous result. A 
small value of A for alkalis (a few opportunities for U-processes) is compensated by 
their quite large value of (T ee )~ l /n, resulting in A ee ’s comparable to those of noble 
metals. One should note that, in some cases, specifically that of potassium and 
copper, the experimental data are strongly sample dependent. This is indicated 
by a wide range of values entered in Table 6. As one goes toward polyvalent metals, 
A ee becomes rather large. Polyvalent metals are often superconductors, and this 
fact seriously curtails the temperature range where the electron-electron interaction 
has a chance to be clearly manifested. Frequently, the data are an admixture of 
electron-electron and electron-phonon interactions, and one must try to separate 
the two by modeling the behavior of the electron-phonon term. We include in 
Table 6 a T 2 term for the resistivity of Nb 3 Sn, one of the highest T, conventional 
superconductors that, close to its transition temperature T c ~ 18 K, displays a 
robust quadratic T-dependence of resistivity that extends to about 40 K. While 
other scattering mechanisms may be at play, it is nevertheless intriguing to consider 
a proposal by Kaveh and Wiser 88 that such a giant T 2 term is not inconsistent with 
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the electron-electron interaction, especially when one takes into account a strong 
enhancement in A ee due to phonon mediation. Surely, in this high-transition-tem- 
perature conventional superconductor, such an enhancement would be particularly 
large even though it must disappear at high-temperatures where the resistivity 
becomes sublinear. There are other features in the transport behavior of these so- 
called A15 compounds that make Kaveh and Wiser argue that the electron-electron 
interaction is a reasonable prospect in spite of the controversy such a proposal 
causes. 

For completeness we also provide in Table 6 the data on two archetypal semi- 
metals, bismuth and graphite. Low and ultralow-temperature measurements show a 
quadratic variation of resistivity with temperature. Whether this is a signature of 
intervalley electron-electron — or, more precisely, electron-hole — scattering 93 or a 
footprint of highly elongated (cigarlike) carrier pockets of electrons on the carrier- 
phonon scattering 116 has been a puzzle and a source of controversy for some time. 
If we take a position that electron-hole scattering is important, it cannot be the 
usual intrapocket variety but rather the interpocket (or intervalley) mechanism. The 
reason is that the carrier pockets in Bi and graphite are very small (carrier densities 
are four orders of magnitude lower than in typical metals) and take up only a small 
fraction of space in the Brillouin zone. The U-processes, the only processes in the 
case of intravalley scattering that can give rise to resistivity, are simply not acces- 
sible in the case of Bi and graphite. On the other hand, interpocket scattering, with 
or without umklapp processes, is in principle a resistive mechanism, just as the case 
of sd scattering in transition metals. In compensated semimetals, i.e., systems with 
an equal number of electrons and holes, even with just N-processes, the interpocket 
scattering is a highly resistive process. 122 This should be evident when one realizes 
that, in response to an electric field, electrons and holes move in opposite direc- 
tions, and thus an interchange of momenta slows down both kinds of carriers. If a 
semimetal is perfectly compensated, the T 2 dependence of resistivity should be 
present at all temperatures. In the case of incomplete compensation, the initial 
T 2 dependence at low temperatures gives way to saturation at high-temperatures. 

It is important to remember that, within the spirit of the relaxation time approx- 
imation, the relaxation time r e - e in Eq. (131) represents an average of the time 
between electron-electron collisions over the entire Fermi surface. However, not all 
collisions are equally effective in hindering the transport, i.e., contributing to the 
resistivity. We mentioned that N-processes cannot give rise to resistivity. One must 
therefore somehow capture this notion of ineffectiveness of some of the collisions, 
and this is accomplished by introducing a parameter A that measures the fraction 
of U-processes among all electron-electron processes. 82 The parameter A plays a 
similar role as the parameter (1 - cos 9 ) — it measures the effectiveness of scattering 
processes in degrading the electrical current. 

In alkali metals normal electron-electron processes are far more frequent than 
the U-processes, and thus A is very small. In noble metals, and especially in poly- 
valent metals, the fraction of U-processes drastically rises and rivals that of N- 
processes. This, of course, reflects the complicated, multisheet Fermi surface pro- 
viding many more opportunities for U-processes. Hitherto one has to be careful 
when calculating coefficients A ee . One cannot take a position that all A ee should be 
roughly the same because all metals have comparable carrier densities and the 
strength of the electron-electron interaction is approximately the same. The struc- 
ture and the shape of the Fermi surface really matters, and through the parameter 
A it has a considerable influence on the transport properties. Moreover, it depends 
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on whether a metal contains isotropic or anisotropic scatterers. The latter situation 
might lead to a considerable enhancement in A ee at low temperatures if N-processes 
dominate because the anisotropic electron scattering makes the N-processes re- 
sistive. 88 

4.3.2. Effect of e-e Processes on Thermal Resistivity 

In the preceding paragraphs we have described a successful approach to determine 
the electron-electron contribution to electrical resistivity. Of course, one might 
hope to use the same strategy to evaluate the coefficient B ee of the electron-electron 
term in the thermal resistivity, Eqs. (136) and (138). Unfortunately, this is not 
possible, at least not by focusing on the low-temperature transport. The difficulty 
arises not just because now all three terms in Eq. (138) are temperature dependent, 

W = B imp T-' + B ee T + f? e -pT 2 , (140) 

but, primarily, because the electron-electron contribution is really small and ther- 
mal transport measurements do not have a prayer of achieving the desired precision 
of a few parts per million. To measure thermal gradients (or temperature differ- 
ences) is a far more challenging task than to measure voltage differences, and the 
precision one can achieve at low temperatures is at best only about 0.1-1%. So, 
experimental attempts to isolate and measure B ee at low temperatures are simply 
futile. 

However, as shown by Laubitz, 123 there is some prospect of accessing and de- 
termining B ee via measurements at high-temperatures. Although it is not easy and 
the utmost care must be exercised to carry the experiments through to their success- 
ful conclusion, there is, nevertheless, a chance. The approach is based on an ob- 
servation that for noble metals the experimentally derived Wiedemann-Franz ratio, 
Tex p = Pexp ( T) / T lUexp ( T ) never quite reaches the theoretical Lorenz number L a = 
(7r 2 /3)(fc B /e) 2 , even at temperatures several times the Debye temperature. The es- 
sential point here is that L exp is lower than L a by a couple of percent as determined 
by the electronic thermal resistivity only; i.e., after any, however small, lattice 
thermal conductivity contribution is subtracted and does not enter consideration. 
Then, because at high-temperatures (well above the Debye temperature) both im- 
purities and phonons scatter electrons purely elastically, these two scattering pro- 
cesses yield the Sonnnerfeld value L a = 2.44x10 8 V 2 /K 2 for the Lorenz number; 
i.e., 

W imp (T) + W e - P (T) = P™r + Pe-P^ . ( 141 ) 

1 L/q 

To isolate a small term IU e _ e against the background of a very large W e _ p , one 
writes, as Laubitz did, the total measured electronic thermal resistivity tl /exp and 
electrical resistivity p exp as 

tyexp p ex P _ WgP(T) pfJ e (T) _ p A ee 
T L a T 2 T L a T 2 ee L 0 ' 1 j 

In principle, Eq. (142) allows for a determination of B ee via high-temperature 
measurements. In reality, because one deals with a very subtle effect, great care 
is needed to account for any small contribution that we otherwise would have 
completely neglected. For instance, we take for granted that the Lorenz ratio is 
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equal to L a = pn/T, while actually, taking into account second-order effects in the 
thermal conductivity [see Eq. (73)], it should be L 0 - S 2 , where S is the thermo- 
power. Under normal circumstances we would completely neglect S 2 because in 
noble metals at ambient temperature it is on the order of 10 11 10 12 V 2 /K 2 (three 
to four orders of magnitude less than L a ) . Likewise, because the experiment is run 
at high but finite temperatures and not at T — > oo, there is a remaining “tail” of 
phonons that scatter electrons inelastically and one should account for it. This is 
accomplished by introducing a term C/T 2 . Thus, in its practical form, Eq. (142) is 
written as 

^exp C ( A \ 

a w^(t) = - {l p _ S 2 )t = Y 2 + [ ■ B <* - ir o J T - ( 143 ) 

Using this approach, Laubitz and his colleagues were able to extract the values of 
B ee for several alkali, noble, and polyvalent metals. Table 7 lists the experimental 
results together with a few calculated values. The unusually large probable errors 
included in the entries in Table 7 should not be surprising, as it is truly a proverbial 
needle in the haystack type of measurement. 

The data in Table 7 indicate good agreement between theoretical and experimen- 
tal values for the alkali metals, whereas calculated values are roughly a factor of 4 
larger than the experimental values for noble metals. That A ee and B ee share a 
common factor {r ee y x /n, 

A ee oc A 
n 



and 



B Pt 



oc • 



(144) 



(the effectiveness parameter A enters only in A ee because for the thermal transport 
all scattering processes hinder heat flow and thus no A is needed) suggests that 



TABLE 7 Values of the coefficient B ee (exp) Given in Units of 10 6 mW 1 Obtained from High- 
Temperature Measurements of the Thermal Conductivity and Electrical Resistivity of Several metals.” 



Metal 


B ee (calc) 


B ee (exp) 


Reference 


Na 


1.1 


1.1±0.6 


34 


K 


3.5 


2.8±0.5 


30 


Rb 


4.9 


3.5±0.5 


41 


Cs 


9.0 


12.4±2 


125 


Cu 


0.22 


0.05T0.02 


123 


Ag 


0.23 


0.05T0.02 


123 


Au 


0.29 


0.09±0.03 


123 


Pb 




~0.5 


126 


A1 




<0.04 


126 



“ Calculated values of B ee (calc) are from Refs. 115 and 121. The entries for Pb and A1 are completely 
unreliable because they are only an estimate and upper bound, respectively. References relate to the 
experimental data. 
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there is a relationship between A ee and B ee . Indeed, Lawrence and Wilkins 82 de- 
rived the following relation: 

Aee = B ee L a g , (145) 

where L a = (7r 2 /3)(fc_B/e) 2 . This relation in principle provides a consistency check 
between A ee and B ee . However, there is a major complication that stems from the 
fact that A ee and B ee are actually temperature dependent, and if a comparison is 
being made it must be done in the same temperature regime — either at very high- 
temperatures or at very low temperatures (but remember, no experimental data are 
available for B ee at low temperatures). The temperature dependence enters because 
of two phenomena : phonon mediation in the electron-electron processes, and the 
influence of anisotropic scattering centers such as dislocations. 

Phonon mediation was introduced by MacDonald 118 and by MacDonald and 
Geldart 124 to account for the fact that the overall electron-electron interaction 
consists of two mutually opposing contributions — the repulsive Coulomb interac- 
tion and the attractive phonon mediated interaction. While the influence of phonon 
mediation differs in different classes of metals (a very small, ~ 1 0%, enhancement in 
noble metals; a factor of 2 decrease in alkali metals; and a huge, more than an 
order of magnitude, enhancement in polyvalent metals), it is effective only at low 
temperatures and does not influence data at high-temperatures, T > 9d • Thus, 
following Kaveh and Wiser 88 , we make a sketch of the behavior of B ee (T ) as a 
function of the reduced temperature, T/9 D (Fig. 6). 

The electrical resistivity term A ee is subject to the same phonon mediation influ- 
ence. However, because it also depends on the collision effectiveness parameter A, 




FIGURE 6 A plot of the temperature dependence of the coefficient B ee [electron-electron contribution to 
the electronic thermal resistivity defined in Eq. (136)] normalized to its high-temperature value B ee { oo). 
The curves represent schematic behavior for a typical alkali metal (K), a typical noble metal (Cu), and a 
polyvalent metal (Al). [After Kaveh and Wiser in Ref. 88] 
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the situation is a bit more complicated. It is well known, and we have indicated so 
in Table 6, that for some metals (potassium and copper in particular) the literature 
data on A ee cover a wide range of values. This has nothing to do with the sensitivity 
or precision of the experimental technique but is due to samples with different 
concentration of defects, especially defects such as dislocations that scatter electrons 
anisotropically. We have pointed out that N-processes are not an impediment to 
electric current. This is strictly true only if electrons do not encounter anisotropic 
scattering centers. Since the samples used in the experiments are often annealed at 
different temperatures and for different durations of time, they contain dislocations 
of usually unspecified, but certainly unequal, concentrations. A problem with dis- 
locations is that they tend to scatter electrons anisotropically. So, how does this 
give rise to an additional temperature dependence of A ee l 
At low temperatures where impurity scattering dominates, one is de facto in the 
anisotropic regime if the sample contains a high density of anisotropic scatterers 
such as dislocations. On the other hand, high-temperature transport is dominated 
by electron-phonon interactions, and at T > d D the scattering is essentially isotro- 
pic. This qualitatively different nature of scattering at low and high-temperatures 
may give rise to an additional (beyond phonon mediation) temperature dependence 
of A ee . Whether this actually happens depends on a particular class of metals and 
on the relative importance of normal and umklapp scattering in each specific situa- 
tion. For instance, in the case of alkali metals, the Fermi surface is essentially 
spherical and is contained well within the first Brillouin zone. In such a case, the 
scattering is dominated by N-processes, the effectiveness parameter A is very small, 
and if only isotropic scattering centers were present the coefficient A ee would vanish 
or be exceedingly small. However, if alkali metals contain high densities of disloca- 
tions, N-processes may turn resistive at low temperatures because scattering of 
electrons on dislocations is anisotropic. This can lead to a dramatic enhancement 
in (A ee ) LT by a few orders of magnitude. On the other hand, if U-processes dom- 
inate the scattering events, the effectiveness parameter A is large (close to unity), 
and whether the electrons do or do not encounter dislocations is not going to have 
an additional substantial influence on the already large A ee . Thus, what gives the 
anisotropic scattering enhancement is the ratio A ams (0)/A lso (0) taken at low tem- 
peratures (marked here as 0 K). The upper bound to the increase in A on account 
of electron-dislocation scattering is in essence given by 1/A^, where A v is the 
umklapp electron-electron value of the effectiveness parameter A. Large A tl r, i.e., 
very frequent U-processes, cannot do much for the enhancement, but small A u, i.e., 
the dominance of N-processes, has a great potential for an enhancement if aniso- 
tropic scattering centers are present. Kaveh and Wiser 127 related the ratio A ows (0)/ 
A !SO ( 0) to the anisotropic p a = p imp + p dis i oc 2* p dis i oc and isotropic 
Po = Pimp + Pdisioc = Pimp limits of scattering. If the preponderance of scattering 
events involves U-processes, there is no anisotropy enhancement. If, however, scat- 
tering is dominated by N-processes, then as the electrons scatter on dislocations 
they make N-processes resistive and an enhancement is realized. Since the para- 
meter A contains information on both umklapp electron-electron and electron- 
dislocation scattering, it follows that if the former one is weak then the electron- 
dislocation scattering can have a great effect on A ee (the case of alkali metals). In 
the opposite case when U-processes dominate, essentially no enhancement arising 
from the electron-dislocation scattering is possible. Kaveh and Wiser provide ex- 
trapolations for various regimes of dislocation density in different classes of metals. 
We conclude this section by sketching in Fig. 7 a possible trend in the temperature 




